
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1991-06 

Implementation of multi-frequency 
modulation with trellis encoding and Viterbi 
decoding using a digital signal processing board 

Wisniewski, John W. 

Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/30967 

This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Callwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkios created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. CaSiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 







AD-A245 998 ' ' 

lllilililliiliilllllllil 

NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



THESIS 


Iinpiementation of Multi-Frequency Modulation 
with Trellis Encoding and Viterbi Decoding 
Using a Digital Signal Processing Board 


by 


John W. Wisniewski 
June 1991 


Itiesls Advisor: 


P.H. Ncnse 


Approved for public release; distribution is unlimited. 


92-03678 









REPORT DOCUMENTATION PAGE 


form Approved 
0MB No 0704-0186 


la REPORT security CLASSIFICATION 




2a. SECURITY CLASSIFICATION AUTHORITY 


2b. DECLASSIFICATION/DOWNGRADING SCHEDULE 


1b RESTRICTIVE MARKINGS 


3. DISTRIBUTION/AVAILABILITY OF REPORT 

Approved for public release; 
distribution is unlimited 


4. PERFORMING ORGANIZATION REPORT NUMBER(S) 


5. MONITORING ORGANIZATION REPORT NUMBER(S) 


6a. NAME OF PERFORMING ORGANIZATION 16b OFFICE SYMBOL I 7a. NAME OF MONITORING ORGANIZATION 

I (If ipplicabte) I 


Naval Postgraduate School 


6c. ADDRESS (City. State, and ZIP Code) 

Monterey, California 93943-5000 


Naval Postgraduate School 


7b. ADDRESS (City, State, and ZIP Code) 

Monterey, California 93943-5000 


8a. NAME Of FUNDING/SPONSORING 
ORGANIZATION 


8b. OFFICE SYMBOL I 9 PROCUREMENT I.MSTRUMENT IDENTIFICATION NUMBER 
(If applicable) | 


Be. ADDRESS (C/ty, State, and ZIP Code) 


10 SOURCE OF FUNDING NUMBERS 


PROGRAM I PROJECT 

ELEMENT NO I NO. 


PROJECT 

TASK 

NO. 

NO 


WORK UNIT 


11 TITLE (Include Security Classification) 

OTLEMENTATION MODULATION WITH TRELLIS ENCODING AND VITERBI DECODING 

USING A DIGITAL SIGNAL PROCESSING BOARD 


12 PERSONAL AUTHOR(S) 

WISNIEWSKI, John William 


13a. TYPE OF REPORT 

Master's Thesis 


13b. TIME COVERED 
FROM_TO 


14. DATE OF REPORT (Year, Month, Day) 

1991 June 


16 SUPPLEMENTARY NOTATION^he views expfessed in this thesis are those of the author and do not 
reflect the official policy or position of the Department of Defense or the U.S. Government. 


17 COSATI CODES 18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 

FIELD I GROUP I SUB-GROUP „_^. „_ _ „ ^ 


Communication; Multi-Frequency Modulation; 
Trellis coding; Viterbi decoding 


19 ABSTRACT {Continue on reverse if necessary and identify by block number) 

Multi-Frequency Modulation Has been the topic of several papers at NPS. In past systems 
the majority of time required for the generation of the MFM signal was due to the software 
routine used to implement the FFT. In this report a Digital Signal Processor was used to 
generate the FFT. The use of Trellis coding and Viterbi decoding on a Digital Signal 
Processor was also investigated. Assembly language programs for three encoder/decoder 
systems were developed. The first uses a 16 QAM signal, the second uses a 2/3 rate 
convolutional encoder and Viterbi decoder and the third uses the CCITT V.32 convolutional 
encoder and Viterbi decoder. 


20 DISTRIBUTION/AVAILABILITY OF ABSTRACT 
iS UNCLASSIFIED/UNLIMITED □ SAME AS RPT 


22a. NAME OF RESPONSIBLE INDIVIDUAL 

MpgjS^P. H. 


DD Form 1473, JUN 86 


21. ABSTRACT SECURITY CLASSIFICATION 
□ OTIC USERS UNCLASSIFIED__ 


22b TELEPHONE (Include Area Code) 22c OFFICE SYMBOL 

‘ EC/Me 




Previous editions are obsolete. 

S/N 0102-LF-014-6603 
i 


SECURITY CLASSIFICATION OF THIS PAGE 

UNCLASSIFIED 









































Approved for public release; distribution is unlimited. 


^ % 


Implementation of Multi-Frequency Modulation 
with Trellis Encoding and Viterbi Decoding 
Using a Digital Signal Processing Board 


by 

John W. Wisniewski 
Lieutenant, United States Navy 
B.S.E.E., Virginia Military Institute 


Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 

from the 


NAVAL POSTGRADUATE SCHOOL 
June 1991 



T. T. Ha, Second Reader 


Michael A. Morgan, Chmrman 
Department of Electrical and Computer Engineering 









0 


ABSTRACT 

Multi-Frequency Modulation has been the topic of several papers at 
NFS. In past systems the majority of time required for the generation of 
the MFM signal was due to the software routine used to implement the FFT. 
In this report a Digital Signal Processor was used to reduce the time 
needed to generate the FFT. The use of Trellis coding and Viterbi 
decoding on a Digital Signal Processor was also investigated. Assembly 
language programs for three encoder/decoder systems were developed. 
The first uses a 16 QAM signal, the second uses a 2/3 rate convolutional 
encoder and Viterbi decoder and the third uses the V.32 convolutional 
encoder and a Viterbi decoder. 
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I. 


INTRODUCTION 


There has long been an interest in frequency division multiplexing as 
a means of combatting impulsive noise, avoiding equalization and making 
fuller use of the bandwidth available. The principle for Multi-Frequency 
Modulation (MFM) was first used over thirty years ago in the Collins 
Kineplex system. Since that time MFM has been used under many names: 

• multiplexed Quadrature Amplitude Modulation (QAM) 

• orthogonal Frequency Division Multiplexing (FDM) 

• dynamically assigned QAM 

• multicarrier modulation 

A parallel or multiplexed data system offers the potential to reduce 
some of the problems of serial systems. In order to increase the data 
rate in a serial system you must utilize higher order modulation or 
decrease the symbol interval, increasing the bandwidth, at the risk of 
degradating the performance of the system. In a parallel data system 
several sequential streams of data are transmitted simultaneously. In a 
classical parallel data system there are N non-overlapping subchannels 
with each data element occupying only a small portion of the available 
bandwidth. 

With the continuing development of Digital Signal Processing 
technology there has been renewed interest in MFM for possible use in: 

• General Switched Telephone Network (GSTN) 

• 60-108 Khz Frequency-division Multiplexed (FDM) group band 
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• Cellular radio 


• High speed data for transmission on the high-rate digital subscriber 

line (HDSL) 

Multi-Frequency Modulation (MFM) can easily be implemented on a 
computer due to the minimal requirements of hardware to construct a 
transmitter or receiver. The primary components of a MFM 
transmitter/receiver are D/A, A/D converters and a method for computing an 
FFT. MFM has been the topic of several papers and projects at NFS using 
various techniques for data acquisition and modulation. One of the 
primary problems with implementing an MFM system in past projects, has 
been the time required to perform an FFT using current software routines. 
With the continuing development of Digital Signal Processing (DSP) chips 
there exist a great number of DSP boards complete with D/A, A/D 
converters, filters, and memory available for industry standard computers. 
This provides the capability to purchase an off the shelf item that has 
all the necessary components to implement an MFM transmitter/receiver. 
Utilizing a DSP board for the computation of the FFT's required for MFM 
provides a tremendous performance gain over the software routines that 
were used in previous projects. Once the programs are downloaded, the 
performance of the DSP board is independent of the performance of the host 
computer. 

General purpose Digital Signal Processor chips, such as the Motorola 
DSP 56001 (which was used in the development of the programs included), 
have an architecture ideally suited to the rapid calculation of FFTs. The 
56001 is a fixed point Digital Signal Processor which has three internal 
memory and address busses which allow the execution of instructions with 
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parallel data moves to the X and Y memory locations. This allows the 
rapid execution of complex computations. There are also 256 locations of 
X and Y Rom that contain Mu-law and A-law expansion tables as well as a 
full four quadrant sine wave table. 

The first objective of this project is to investigate the use of a DSP 
for implementing an MFM system and reduce the time required to generate a 
MFM signal. The second objective was to investigate the use of trellis 
coding and Viterbi decoding using a Digital Signal Processor. 
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II. THEORY 


A. MULTI-FREQUENCY MODULATION 


A Multi-Frequency Modulated signal has a packet structure and is 
comprised of time and frequency slots.(Figure 1) 



Figure 1 MFM Signal Packet (from Ref. 1: p.3) 


The following definitions will be used for discussing MFM: [Ref. l:pp 5-6] 
• K: Number of MFM tones 


• T: Packet length in seconds 






• aT: Baud length in seconds 

• k^: Baud length in number of samples 

• L: Number of bauds per packet 

• At: Time between samples in seconds 

• f/=l/At: Sampling or clock frequency for D/A and A/D conversion in 
H2 

• Af=2/AT: Frequency spacing (minimum) between MFM tones 

• Symbol set. Phase of the k^** tone in the l^** baud 

• Amplitude of the k^** tone in the l^** baud 

Each MFM packet consists of L bauds of K tones. The packet 

construction for the 1^** baud is given by: [Ref 1: pp 6-7] 

Xj(u)-^;cj*(u) (1) 

where, the analog representation of each tone during the l^*" baud is given 
by: 

The time ,u , is referenced from the beginning of the baud. The time at 

the beginning of the baud is defined as tp and the time at any given point 

in the packet is t=tQ+lAT+u. 

By sampling (1) and (2) at intervals At^l/fj^ a discrete time sampled 
version for the 1^ baud can be found: 

(3) 


5 


where the discrete time version of each tone during the baud is given 
by: 

Xjif in) ~ Ai^coB (2nkn /; Osnskjf-1 (4) 

The value, n, is the discrete time referenced to the beginning of the 
baud. From the Sampling Theorem the maximum frequency must be less than 
f^/2. There are a maximum of k^/2 tones available spaced at intervals of 
Af between dc and f^/2-Af. The value of k is in the range from 0 and kjj/2, 
where k is the harmonic number. 

By taking the Discrete Fourier Transform (DFT) of (3): 

x,ik)-'£x,,ik) (5) 

The discrete time signal (3) can be generated by the k^ point IDFT of (5): 

x^{n)-IDFT (6) 

It is easily seen that the l^*" baud is generated by taking the IDFT of a 

complex valued array of length k^j. The values in the array, Xj(n) are the 
discrete time samples of the analog transmit signal. The generation of 
the 1^** baud is completed by clocking the k^^ samples out at fj, samples/sec. 
The entire packet is completed by an L fold repetition of the procedure. 

B. TRELLIS CODED MODULATION 

Trellis coded modulation evolved as a combination of coding and 
modulation techniques for digital transmission over band limited channels. 
The primary advantage for using trellis coded modulation over other coding 
schemes is that significant coding gains can be achieved without 
compromising bandwidth efficiency [Ref 2:pp 5-12]. Trellis coded 
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modulation schemes employ redundant nonbinary modulation in combination 
with a finite-state encoder to govern the selection of modulation signals 
and to generate coded signal sequences. A simple four state scheme can 
improve the reception of a digital transmission by as much 3db in additive 
white gaussian noise. If more complex coding schemes are used gains of 
6db or more may be achieved. In order to achieve the potential gains of 
a trellis coded signal a soft decision decoder, such as a Viterbi decoder, 
must be utilized in the receiver. These gains can also be achieved 
without reducing the information rate or increasing the bandwidth as 
required by conventional error correction schemes. 

For low to medium data rates (s 4800 bits/s), signaling methods that 
use independent symbol-by-symbol transmission are adequate to provide 
acceptably low error rates over voice grade circuits. When data rates 
increase in speed (2: 9600 bits/s) the same is not true, QAM and optimal 
signal sets fail to provide acceptably low error rates.[Ref 3:pp 648-649] 

There are many factors which comprise linear and nonlinear distortion 
on a voiceband channel including phase jitter, frequency offset and 
additive noise. It was shown in that, relative to the uncoded system, 
trellis codes with four and eight states provide marked improvement in 
performance with respect to additive noise, second and third harmonic 
distortion, phase jitter, impulse noise and other channel imparements.[Ref 
3:pp 649-651] 

For an uncoded system the binary data transmit rate is equal to m/T 
bits/sec, where m is the number of bits/symbol and 1/T the symbol rate. 
In a conventional QAM system (uncoded) there are 2" discrete symbol points 
(amplitude and/or phase levels) with each successive symbol transmitted 
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independently. The error performance depends on the minimum distance 
between the signal points (the larger the distance, the lower the error 
rate). The minimum distance at the transmitter is limited by the average 
power allowed on the circuit and the choice of the signal points. Even 
with the use of in-phase and quadrature components the minimum distance 
(d^j^g) between points decreases with the increase in the number of symbols, 
m, and constant average power. This results in a degradation at higher 
transmission rates assuming that there is a constant symbol rate. 

The objective of coding is to increase the effective minimum distance 
between the signals without increasing the average power. One method of 
accomplishing this is with the use of convolutional encoder. In a 
convolutional encoder there are m current bits, and v past bits used to 
develop the codeword. The vpast bits define the state of the encoder and 
are operated on to produce a+j bits, the rate of the code is described as 
m/m'tj. The m+j bits require 2“*^ discrete channel symbols. Using a 
convolutional encoder the minimum distance between symbols is no longer 
the measure of error performance, performance is now a measure of the 
minimum distance between the allowed transition of symbols from one state 
to another.[Ref 3:p 649] A convolutional encoder using a shift register 
to provide two past bits for a 2/3 rate convolutional code is shown in 
Figure 2. The allowed transition between states is shown in Figure 3. A 
convolutional code can be described as an (n,k,m) code, where n is the 
number of encoded bits, k the number of information bits and m the number 
of past bits used for encoding. 
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Figure 2 2/3 Rate Convolutional 
Encoder (from Ref 3:pp 651) 

The input and output relationship are depicted by the branches on the 
trellis diagram. The upper branch depicts the transition with Xq set to 
"0”, the lower branch depicts the transition with Xq set to "1". Each 
node of the trellis diagram represents one of four states created by the 
past bits Sj and S2. 



Figure 3 Trellis for 2/3 rate 
convolutional code.(from Ref 3:pp 651) 


In order to achieve the optimum decoding gains from the use of 
convolutional encoders the decoder must use a trace back routine to find 







the most probable path through the trellis. The rules for bit to symbol 
mapping for coded systems:[Ref 3:pp 650] 

1. All parallel transitions in the trellis structure receive maximum 
possible Euclidean distance in the signal constellation. 

2. All transitions diverging from a merging into a trellis state receive 
maximum possible Euclidean distance. 

C. VITERBI DECODER 

The Viterbi decoder utilizes a maximum likelihood sequence estimation 
method to decode the incoming data stream. A predetermined measure is 
used to determine the symbol sequence which is closest to the received 
symbol sequence. At any time, k, the shortest path, called the survivor, 
entering each state (node) of the trellis is retained. To proceed to time 
k-fl, all time k survivors are extended by computing the metrics (lengths) 
of the extended path segments based on the calculated branch metrics, 
dependent on the branch symbols in the trellis and the value of the 
received sample.[Ref 3:pp 648-649] An example of a trellis and a path 
through the trellis is shown in Figure 4. The metrics of the remaining 
paths are computed and the shortest path retained. The shortest length 
into each state, the (k4^1) survivor is retained. The number of survivors 
never exceeds the number of states in the trellis. 

The basic operations required in a soft decision Viterbi decoder 
are:[Ref 3:p 651] 

* computation of branch metrics, for p additive white gaussian noise 
channel, are proportional to (r-x^)^ where r is the received sample 
and x^ is the noise free symbol associated with the message. 

• addition of branch metrics and survivors to determine the survivors 
for each state. 
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Figure 4 Trellis diagram for a (3,1,2) code, (from Ref 4:p 
316) 

• comparison along the extended path metrics to determine the survivor 
for each state. 

It was shown that the Viterbi algorithm was equivalent to a dynamic 
programming solution to the problem of finding the shortest path through 
a weighted graph [Ref 3:pp 317-321]. The decoder must produce an estimate 
P of the codeword v based on the received sequence r. For an information 
sequence of length L, the trellis must contain L-^m-fl levels or time units 
to decode the sequence. An (n,k,m) code has an information sequence of 
length Kl, and is encoded into a codeword of length Nsn(L-«-m). A maximum 
likelihood decoder (MLD) for a discrete memor/less channel (DMC) chooses 





as the codeword v which maximizes the log-likelihood function log 
P(r|v). Since for a DMC: 

tt-i 

I -JJ P(r J Vj) (7) 

it follows that log-likelihood function is formed by the summation of the 
branch metrics: 

L*m~x N-l 

logP(r I v) - ^ logP(r^ I v^) logP(rj( I v_() (8) 

where P(r.|vj) is a channel transition probability. This is a minimum 
error probability decoding rule when all code words are equally likely. 
The log-likelihood function log P(r|v) are called the metric of path v and 
is denoted M(r|v). A partial path metric is formed by summing up the 
partial path metrics for J branches and is expressed by: 

(Q\ 

M{ [r I v] j) I v^) 

The final survivor in the Viterbi algorithm is the maximum likelihood 
path; 

Af(r I V) ^Af(r I V) , all V ^ ^ (10) 

The path is then traced back to determine the symbol transmitted.[Ref. 
.4:pp 316-318] 


12 



III. SYSTEM DEVELOPMENT 

All the programs used in this thesis were run on an Ariel PC-56 DSP 
board. The PC-56 DSP board is an eight bit card that can be used in any 
industry standard computer. The primary components of the PC-56 are a 20 
Mhz Motorola DSP 56001, 16k of external memory, a 24 bit bidirectional 
interrupt-driven port with a header for external I/O, a TLC32040 14 bit 
analog interface chip with built in input and output filters gain section. 
An Ariel PC-56D was also used for comparison, it contains a 27 Mhz 
Motorola DSP 56001, 64k of external memory and a NEXT compatible DB-15 
port. For actual implementation of this system a DSP board with dual A/D 
and D/A converters is required. 

The structure of the Motorola DSP 56001 makes it ideally suited for 
high speed communications. In the DSP chip there are 256 locations of 24 
bit X and y memory that occupy the lowest 256 locations of the DSP address 
space. The locations from 256-511 are allocated for the on-chip ROM. The 
onboard ROM contains Mu-law and A-law expansion tables as well as a full 
four quadrant sine table. There are 512 locations of 24 bit high speed 
program RAM (PRAM) on the chip.[Ref 5] 

A feature that makes the DSP 56001 desirable for use in this system 
is that it offers several addressing modes. It implements three types of 
arithmetic for addressing, linear, modulo and reverse carry. For each of 
the address registers R0-R7 there is an offset register N0~N7 and a 
modifier register M0-M7. The offset register contains the values to 
increment and decrement the address register. The modifier register 
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defines the type of address arithmetic to be used. For modulo arithmetic 
the contents of the modifier register Mn specify the base modulus.[Ref 
6;pp 5-2,5-4] Tha DSP 56001 also provides three different addressing 
modes, register direct, address register indirect, and post or pre 
increment/decrement. 

The basic configuration of an HFM transmitter and receiver are shovm 
in Figure 5 and Figure 6. For the purpose of this project three types of 
modulators were used in the investigation for use on a DSP board. In 
previous systems implemented at NPS [Ref 1] and [Ref 7] the complex 
conjugate of the input data was loaded into the image frequencies of the 
IFFT. This resulted in only real data being generated by the IFFT. For 
this project complex data was loaded into all available locations to 
generate both real and complex data. The signal ma^ be modulated using 
any type of modulation scheme. For the purpose of this thesis the three 
different modulators used a 16 QAM signal, a 2/3 rate convolutional 
encoder with an eight PSK signal and a 4/5 rate code with a 32 QAM signal 
using the CCITT modem standard V.32 convolutional encoder. The source 
code for the V.32 encoder and Viterbi decoder was included in an example 
manual from Motorola [Ref 8:pp A1-C2] for use on their simulator and was 
modified for use on the Ariel board to generate a Multi-Frequency 
Modulated signal. 
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Figure 5 MFM Transmitter 



Figure 6 MFM Receiver 
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A. ENCODERS 


Of the 256 bins available in the IFFT array, only 201 message symbols 
are carried in each baud. The other 55 bins are loaded with 2 eros to 
allow for filtering of the signal. Once the 256 bins are filled, the IFFT 
subroutine is called to generate the values for the MFM signal. At this 
point the values would be clocked out through the D/A converters. 

1. 16 QAM Encoder 

At initialization of the system the message to be transmitted is 
downloaded from the host computer as well as the programs for encoding 
(QAK16EN) the signal and the data for the look-up table (16QAMRE.DAT and 
16QAMIM.DAT). The data is read in using four bit increments. The value 
read in is moved into an offset register and used to determine the 
coordinates of the points on the constellation in Figure 7. The real 
values of the constellation are stored in the x. memory and the 
corresponding imaginary values are stored in the y memory. Once the 
constellation values are determined the real and complex values are stored 
in memory. Once one hundred samples have been read into memory the 
routine installs 55 zeros at the center of the array then reads the 
remainder of the 256 symbols. Once all 256 samples are stored the IFFT 
routine is called to generate the time domain signal. 

2. 2/3 Sate Convolutional Encoder 

At initialization of the system the message to be sent as well as the 
encoder (23ENC0D) and the data files (23REAL.DAT and 23IHAG.DAT) are 

I 

loaded into memory. In the 2/3 rate encoder two bits are read in from the 
input message and using the convolutional encoder of Figure 2 a third bit 
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Figure 7 CONSTELLATION FOR 16QAM SIGNAL (from 
Ref.8:p 2-2) 


is generated. The eight PSK constellation in Figure 8 is used to transmit 
the message. A look up table is used to find the points on the 
constellation. After the points on the constellation are found, the real 



Figure 8 2/3 Rate Code 

Constellation 
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and imaginary values are loaded into an array. After the first 100 
samples are loaded into the array, 55 zeros are loaded into the array to 
allow for filtering and the remainder of symbols are read. Once all 256 
values are loaded the IFFT is performed and the values are stored ready to 
be clocked out. 

3. V.32 Convolutional Encoder 

At initialization the message, encoder program (MFMENCOD) and data 
files (QANREAL.DAT and QAM1MAG.DAT) are downloaded to the to the DSP 
board. For the V.32 encoder after the four bit ssnnbol is read, the 
convolutional encoder Figure 9 is used to generate a fifth redundant bit. 



Ref. 8:p 2-4) 


The 32 QAM constellation in Figure 10 is used for the five bit signal. A 
look up table is used to determine the real and imaginary values for the 
corresponding point on the constellation. As in the 16 QAM encoder and 
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the 2/3 rate code encoder, 201 encoded symbols and 55 zeros are loaded 
into the IFFT array and the IFFT is performed. 

B. GENERATING THE MFH SIGNAL 

When all 256 bins have been loaded the program calls the routine to 
perform the IFFT. The routine used to perform the IFFT was generated by 
using the relationship x(n)s(FFT[X*(k)]*/N. The program to generate 
the FFT was included with the Ariel DSP packaged In order to speed up 
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Figure 10 32 QAM Constellation 
(from Ref. 8:p 2-3) 


Unless otherwise specified, all software (C) COPYRIGHT 1989 by Ariel 
Corporation, Highland Park, NJ. All rights reserved. You may use the 
software provided by Ariel Corporation on any one computer; copy the 
object code into my machine readable form for your use. You may modify 
the software provided by Ariel and/or merge or incorporate it into any 
general use software program of your development except for a program of 
similar nature to the Ariel product. You may freely reproduce any such 
program of your development; however the merged or incorporated part of 
the Ariel provided software will continue to be subject to all other 
provisions of this agreement. 
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the routine as much as possible the values used in the look up table were 
conjugated and scaled by 256 to save from doing them in the modulator. By 
using the prescaled and conjugated values in the encoder constellations 
the IFFT routines require no more time than the FFT routine used in the 
decoder. 

After computing the IFFT the real and complex values are stored in 
memory. By utilizing a DSP board with dual D/A converters the real and 
imaginary values can be clocked out at the sampling frequency f^^. The 
design of the system may be altered to either allow transmission of each 
256 symbol packet or store the entire message and transmit at the end of 
the message. 

C. DECODERS 

After using a quadrature receiver to recover the inphase and 
quadrature components, the signals are sampled at the same frequency, f^^, 
as used in the encoder. These sampled values are used to determine the 
transmitted symbols. 

1. 16 QAH DECODER 

At initialization the boundary data file B0UN’16.D and the decoder 
program QAM16DEC must be loaded onto the DSP board. To decode the 16 QAM 
signal the constellation was partitioned as shown in Figure 11. The 
magnitude of each received signal value is compared to the boundaries in 
the first quadrant. Once a bounded area is found the actual values are 
then used to determine the correct quadrant. The data file B0UN16.D is 
indexed and the corresponding symbol is read out. 
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Figure 11 PARTITIONED 16 QAM CONSTELLATION 


2. 2/3 Rate Code Decoder 

The use of a convolutional encoder requires the use of a soft decoder 
such as a Viterbi Decoder to recover the transmitted symbol. The boundary 
data file BOUN23.D and the Viterbi decoder program KFM23DEC are first 
loaded into memory and run. After computing the FFT to recover the 
transmitted symbols, a received point is read out of the array and the 
quadrant of the constellation (Figure 8) that it lies in is found. The 
data file B0UN23.D holds the four closest points, one in each state, to 
the quadrant that the received point lies in. A boundary file is used to 
help speed up the decoder routine, instead of computing the distance to 
all points in the constellation only the point closest in each state is 
used for computing the euclidean distance to the received point. Using the 





received point and the four values found in BOUN23.D the euclidean 
distance to each point is computed. The accumulated distance table is 
updated with the distance to each state. The smallest accumulated 
distance is found and the trellis is traced back 16 time periods to find 
the most likely transmitted point. After finding the most likely point it 
must be decoded. The bit Y2 is discarded and using the relation YlsXl$Sl 
the value of XI is recovered by using past state information and Xl=yi$Sl. 
The value of YO is output as XO. The values are stored in memory. 

3. V.32 DECODER 

The use of a convolutional encoder in the V.32 system requires that 
a Viterbi Decoder must also be used to recover the transmitted symbol. 
The boundary file BOUND.D and the program file MFMDECOD are first leaded 
into memory. The constellation must first be partitioned in such a way 
that the partition equally divides the distance between the four symbols 
in the same state. The partition used for the state 110 is shown in 
Figure 12. By using this method the partitions for all eight states may be 
imposed on the constellation resulting in 52 separate partitions as shown 
in Figure 13. For each partition there are eight points, one from each 
state, that is closest to the boundary. These eight points for each of 
the 52 boundaries are stored in the data file BOUND.D that is loaded into 
tlie processor at initialization. Once a bounded area has been found, the 
eight points are used to determine the euclidean distance from the 
received point to the points read from BOUND.D. The distances calculated 
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P 3-7) 



Constellation (from Ref. 8:p 
3-7) 


are used to update the accumulated distance to each state. The state with 
the smallest accumulated distance is found and the trellis is traced back 
16 time periods to find the input state at the end of the most likely 
path. The closest point in that state to the input point at that time 
period is found and output. The most significant bit is masked off. The 
differential encoding used on the two most significant bits must be 
decoded by: 

Once the two most significant bits have been decoded the symbol is output 
to memory. 
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IV. SYSTEM OPERATION 


The Digital Signal Processing boards used to run the programs that 
were developed did not have dual D/A, A/D converters required to fully 
implement an HFM system. The programs as listed, do not fully implement 
a MFM system. There are no subroutines included for the clocking out of 
samples through D/A converters in the encoders and no sampling routines in 
the decoders. The encoding programs read in the data to be transmitted 
from memory, encode the data and store the values created by the IFFT in 
memory. The same locations are used to store all values generated, it is 
assumed that the values will be clocked out once available freeing up 
these locations for the next baud. The values would be clocked out through 
the D/A converters to generate a signal. In the decoder programs, the 
data is read out of locations were the sampled values would be stored. 
The same locations are used for all sampled values, it is assumed each 
baud will be decoded prior to receiving the next baud. If a convolutional 
encoder is used, a decision must be made as to whether the next baud 
should be delayed to allow decoding of the previous baud, or transmit the 
baud when available and store the values in the decoder for off line 
decoding. 

Once the modulation technique is chosen the appropriate constellation 
data files must be loaded for the look up tables used in the decoder. 
These files are located in Appendices D, H, L. Once the appropriate data 
files are loaded the encoder program and message file are loaded into 
memory and the encoder is run. The programs developed were run out of 
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Ariel's BUG-56 monitor-debugger, they can also be adapted to run as 
subroutines in a Microsoft C program. The monitor-debugger allowed easy 
access to all the chips registers and functions making running and 
debugging the programs simpler. 

To use the decoder the appropriate boundary values must be loaded. 
Once the boundary values are loaded the decoder program is loaded and run. 
The output is stored in memory. 

All of the programs listed are written in assembly language for the 
DSP 56001 and must be compiled into a loadable file using the Motorola 
assembler. All of the data files are properly formatted to be read into 
the correct memory locations. If the programs are run out of the BUG-56 
monitor-debugger, a macro can be created to load all the files needed for 
an encoder or decoder using one instruction. 
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V. CONCLUSIONS 


The time required for the encoder/decoder for a 256 sample baud is 

shown in Table I. The Turbo Pascal routine used in past projects 

Table I Sample Times for Programs Using 256 
Sample Baud 


TIKE(ms) 

PROGRAKS 20 Khz 56001 27 Khz 56001 


4096 pt complex fft .023 sec .017 sec 

256 pt complex fft 1 .7 


MFM PROGRAMS (for 256 sample baud) 


16 QAK Encoder 

2.7 

2.0 

16 QAM Decoder 

2.7 

2.0 

2/3 convolutional 
encoder 

3.2 

2.4 

2/3 rate code 

Viterbi decoder 

10.8 

8.0 

V.32 encoder 
(w/o diff enc) 

3.6 

2.7 

V.32 Viterbi decoder 
(w/o diff enc) 

15.8 

11.7 

V.32 encoder 

3.7 

2.8 


V.32 Viterbi decoder 


16.9 


12.5 






required 22 seconds to perform a 4096 point complex FFT on an 8 Mhz AT 
computer with a coprocessor.[Ref• IB] It is readily apparent that the 
Digital Signal Processor has a tremendous performance advantage over 
normal software routine. Using Matlab benchmarks a Sun Sparc station 
required .34 seconds and a 20 Mhz 386 with a coprocessor required 1.16 
seconds for a 4096 point complex FFT. 

Multi-Frequency Modulation can be easily implemented on a Digital 
Signal Processor Board. The programs for implementing the 16 QAM system 
required the same amount of time to encode and to decode. There is little 
penalty for using a trellis encoder but the Viterbi decoder takes over 
three times longer than the encoder for the 2/3 rate code and almost four 
times as long for the V.32 decoder. 

The use of DSP boards in implementing MFM shows much improvement over 
existing techniques. Suggestions for future research include running the 
MFM programs on DSP boards with the required dual D/A, A/D converters. A 
determination must be made as to whether the decoder should store the 
entire message or delay transmission of a baud until the previous baud is 
decoded. Another alternative is to use more than one processor, one to 
encode/decode and another to clock out/sample the data. 

The versatility and speed of Digital Signal Processors make them 
ideally suited for use for a Multi-Frequency Modulated System. They offer 
much greater performance than is currently available using software 
routines and are can easily perform other functions by simply loading a 
new program. The continuing development of more DSP compatible products 
show that there will be greater uses of DSP technology in the future. 
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APPENDIX A. 


FFT 


; This program originally available on the Motorola DSP bulletin board. 
; It is provided under a DISCLAIMER OF WARRANTY available from 
; Motorola DSP Operation, 6501 Vm. Cannon Drive V., Austin, Tx., 78735. 

f 

; Radix 2, In-Place, Decimatlon-In-Time FFT (fast). 

; Last Update 18 Aug 88 Version 1.0 
• 

fft macro points,data,odata.coef 

fft ident 1,0 

; Radix 2 Decimation in Time In-Place Fast Fourier Transform Routine 

; Complex input and output data 

; Real data in X memory 

; Imaginary data in Y memory 

; Normally ordered input data 

; Normally ordered output data 

; Coefficient lookup table 

; -Cosine values in X memory 

; -Sine values in Y memory 

; Macro Call - fft points,data,odata,coef 

; points number of points (16-32768, power of 2) 

; data start of data buffer 

; odata start of output data buffer 

; coef start of sine/cosine table 

; Alters Data ALU Registers 


xl 

xO 

yi 

yO 

a2 

al 

aO 

a 

b2 

bl 

bO 

b 


Alters Address Registers 


rO 

nO 

mO 

rl 

nl 

ml 


n2 


r4 

n4 

m4 

r5 

n5 

m5 

r6 

n6 

m6 
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Alters Program Control Registers 
pc sr 

Uses 6 locations on System Stack 
Latest Revision - 18 Aug-88 


move //data,rO ;Initialize input pointer 

move //points/4,nO initialize input and output pointers offset 

move n0,n4 ; 

move n0,n6 initialize coefficient offset 

move //points-l,mO initialize address modifiers 

move mO,ml ;for modulo addressing 

move m0,m4 
move m0,m5 

Do first and second Radix 2 FFT passes, combined as 4-point butterflies 


move 


x:(rO)+nO,xO 



tfr 

xO,a 

x:(rO)+nO,yl 



do 

nO,_i 

twopass 



tfr 

yl,b 

x;(rO)+nO,yO 



add 

yO,a 

x: (rO),xl 


;ar+cr 

add 

xl,b 

r0,r4 


;br+dr 

add 

a,b 

(rO)+nO 


;ar'=(ar+cr)+(br+dr) 

subl 

b,a 

b,x:(rO)+nO 


;br'=(ar+cr)-{br+dr) 

tfr 

xO,a 

,x0 

y:(rO),b 


sub 

yO,a 


y:{r4)+n4,y0 

1 

O 

sub 

yo,b 

xO,x:(rO) 


;bi-di 

add 

a,b 


y; (rO)+nO,xO 

;cr’s:(ar-cr) + (bi-di) 

subl 

b,a 

b,x:(rO) 


;dr'=(ar-cr)-{bi-di) 

tfr 

xO,a 

,x0 

y:{r4),b 


add 

yo,a 


y: (rO)+nO,yO 

;bi+di 

add 

yO,b 

xO,x:(rO)+nO 


;ai+ci 

add 

b,a 


y; (rO) + ,xO 

;ai'=(ai+ci)+{bi+di) 

subl 

a,b 


a,y:(r4)+n4 

;bi’=(ai+ci)-(bi+di) 

tfr 

xO,a 


b,y:(r4)+n4 


sub 

yo,a 

xl,b 


;ai-ci 

sub 

yl,b 

x;(rO)+nO,xO 


;dr-br 

add 

a,b 

x;{rO)+nO,yl 


;ci'*(ai-ci)+(dr-br) 

subl 

b,a 


b,y:(r4)+n4 

;di'=(ai-ci)-(dr-br) 

tfr 

xO,a 


,y:(r4)+ 



twopass 

Perform all next FFT passes except last pass with triple nested DO loop 

move //points/8,nl initialize butterflies per group 

move //4,n2 initialize groups per pass 

move //-I,m2 iinear addressing for r2 
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move ^/0,m6 


;initialize C address modifier for 
;reverse carry (bit-reversed) addressing 


do (/@cvi (@log( points) /@log( 2) -2.5) ,_endjpass 

{example: 7 passes 

for 1024 pt. FFT 
move (/data.rO 


{initialize A input 

pointer 

move r0,rl 
move nl,r2 
move r0,r4 


{initialize A 

output pointer 
move {rl)+nl 


{initialize 6 input 

pointer 

move rl,r5 


{initialize B 

output pointer 

move <(coef,r6 


{initialize C input 

pointer 

lua (r2)+,n0 


{initialize pointer 

offsets 

move n0,n4 
move n0,n5 
move (r2)- 


{butterfly loop 

count 

move x:(rl),xl 

y:{r6),y0 

{lookup -sine and 

-cosine values 

move x:(r6)+n6,x0 

y:(r0),b 

{update C pointer, 

preload data 

mac xl,y0,b 

y:{rl)+,yl 


macr -x0,yl,b 

y:(r0),a 


do n2,_end_grp 
do r2,_end_bfy 
subl b,a x:(rO),b 

,y:(r4) 

{Radix 2 DIT 

butterfly kernel 

mac -xl,x0,b x:(r0)-f,a 

,y:(r5) 


macr -yl,y0,b x:(rl),xl 
subl b,a b,x:(r4)-f 

y:(r0),b 


mac xl,y0,b 

y:(rl)+,yl 


macr -x0,yl,b a,x:(r5)4 

y:{r0),a 


_end_bfy 

move (rl)+nl 

subl b,a x:(rO),b 

,y:{r4) 


mac -xl,x0,b x: (rO)-<-nO,a 

,y:(r5) 


macr -yl,y0,b x:(rl),xl 

y:(r6),y0 


subl b,a b,x:(r4)-«-n4 

y:(r0)»b 


mac xl,y0,b x: (r6)-^n6,x0 

y:{rl)+,yl 


macr -x0,yl,b a,x:(r5)4n5 

y:{r0),a 


_end_grp 

move nl,bl 

Isr b n2,al {divide 

butterflies per group by two 

Isl a bl,nl {multiply groups per pass 

by two 
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move al»n2 
_end_pass 

I 

; Do last FFT pass 

f 

move #2,n0 ;initialize pointer offsets 

move nO,nl 

move #points/4,n4 ;output pointer A offset 

move n4,n5 ;output pointer B offset 

move )'/data,rO ;initialize A input pointer 

move j/odata,r4 ;initialize A output pointer 

move r4,r2 ;save A output pointer 

lua (rO)-f,rl ;initialize B input pointer 

lua (r2)+n2,r5 ;initiali 2 e B output pointer 

move //0,m4 ;bit-reversed addressing for output ptr. A 

move m4,m5 ;bit-reversed addressing for output ptr. B 

move f/coef,r6 ;initialize C input pointer 

move {r5)-n5 ;predecrement output pointer 

move x:(rl),xl y:(r6),y0 

move x:(r5),a y:(rO),b 

do n2,_lastpass 

mac xl,yO,b x:(r6)+n6,x0 y:(rl)+nl,yl ;Radix 2 DIT butterfly 
kernel 

macr -xO,yl,b a,x:(r5)+n5 y:(rO),a ;with one butterfly per 

group 

subl b,a x:(rO),b b,y:{r4) 

mac -xl,xO,b x:(rO)+nO,a a,y:{r5) 

macr -yl,yO,b x:(rl),xl y:(r6),y0 

subl b,a b,x:(r4)+n4 y;(rO),b 

_lastpass 

move a,x:(r5)+n5 

endm 




APPENDIX B. 


IFFT 


This program uses the same routine used in the program FFT. The IFFT is 
calculated using the relationship IFFT=(FFT[X(k)]/N). It is assumed that 
the values used in the constellation are conjugated and prescaled by N. 
If the values to be used are not conjugated a routine must be added to 
the beginning of this routine to conjugate the values. The division by 
N can also be accomplised by using left shifts, but will slow down the 
program. 

Radix 2, In-Place, Decimation-In-Time IFFT (fast). 


ifft macro points,data,odata,coef 

ifft ident 1,0 

; Radix 2 Decimation in Time In-Place Inverse Fast Fourier Transform 
; Routine 

; Complex input and output data 

; Real data in X memory 

; Conjugated and prescaled imaginary data in Y memory 

; Normally ordered input data 

; Normally ordered output data 

; Coefficient lookup table 

; -Cosine values in X memory 

; -Sine values in Y memory 


Macro Call - 

ifft 

points 

,data,odata,coef 

points 


number of points (16-32768, power of 2) 

data 


start of 

data buffer 

odata 


start of 

output data buffer 

coef 


start of sine/cosine table 

Alters Data ALU Registers 


xl 

xO 

yl 

yo 

a2 

al 

aO 

a 

b2 

bl 

bO 

b 

Alters Address Registers 


rO 

nO 

mO 


rl 

nl 

ml 



n2 



r4 

n4 

m4 


r5 

nS 

m5 
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r6 n6 m6 


Alters Program Control Registers 
pc sr 

Uses 6 locations on System Stack 


Latest Revision - 18 Aug-88 


move //data.rO 
move //points/^.nO 
move n0,n4 
move n0,n6 
move //points-l,mO 
move mO,ml 
move m0tm4 
move m0,m5 


;initialize input pointer 
;initialize input and output pointers offset 

;initialize coefficient offset 
;initialize address modifiers 
;for modulo addressing 


Do first and second Radix 2 IFFT passes, combined a'; 4-point butterflies 


move 


x: (rO)+nO,xO 



tfr 

xO,a 

x: (rO)+nO,yl 



do 

nO,_l 

twopass 



tfr 

ylfb 

x:{rO)+nO,yO 



add 

yO,a 

x:(rO),xl 


;ar+cr 

add 

xl,b 

r0,r4 


;br+dr 

add 

a,b 

(rO)+nO 


;ar'=(ar+cr)+(br+dr) 

subl 

b,a 

b,x:(rO)+nO 


;br' = (ar+(:r)-(br+dr) 

tfr 

xO,a 

a,xO 

y:(rO),b 


sub 

yO,a 


y; (r4)4n4,y0 

;ar-cr 

sub 

yO,b 

xO,x:(rO) 


;bi-di 

add 

a,b 


y: (rO)+nO,xO 

;cr'={ar-cr)+(bi-di) 

subl 

b,a 

b,x:(rO) 


;dr'=(ar-cr)-(bi-di) 

tfr 

xO,a 

a,xO 

y:{r4),b 


add 

yO,a 


y: (rO)+nO,yO 

;bi+di 

add 

yo.b 

xO,x:(rO)^nO 


;ai-fci 

add 

b,a 


y:(rO)+,xO 

;ai'=(ai+ci)+(bi+di) 

subl 

a,b 


Styj (r'i)+n4 

;bi'=(ai+ci)-(bi+di) 

tfr 

xO,a 


b,y:(r4)+n4 


sub 

yO,a 

xl,b 


;ai-ci 

sub 

yl,b 

x:(r0)4n0,x0 


;dr-br 

add 

a,b 

x:{rO)+nO,yl 


;ci'a(ai-ci)+(dr-br) 

subl 

b,a 


b,y:(r4)+n4 

;di’e{ai-ci)-(dr-br) 

tfr 

xO,a 


a,y:{r4)+ 



_twopass 

» 

; Perform all next IFFT passes except last pass with triple nested DO loop 
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move 

//points/8, nl 

;initialize butterflies per group 

move 

//4,n2 

;initialize groups per pass 

move 

//-I, m2 

;linear addressing for r2 

move 

//0,m6 

;initialize C address modifier for 



;reverse carry (bit-reversed) addressing 

do 

//0cvi(@log(points)/@log(2)-2.5),_endjpass ;example: 7 passes 

for 1024 pt. IFFT 


move 

i/data, rO 

;initialize A input 

pointer 



move 

r0,rl 


move 

nl,r2 


move 

r0,r4 

;initialize A output 

pointer 



move 

(rl)+nl 

;initialize B input 

pointer 



move rl,r5 

;initialize B output 

pointer 



move 

//coef ,r6 

;initialize C input 

pointer 



lua 

(r2)+,n0 

;initialize pointer 

offsets 



move 

n0,n4 


move 

n0,n5 


move 

(r2)- 

;butterfly loop count 

move 

x:(rl),xl 

y;(r6),y0 ;iookup -sine and -cosine 

values 



move 

x; (r6)+n6,: 

xO y:(r0),b ;update C pointer, preload 

data 



mac 

xl,y0,b 

y:(rl)+,yl 

macr 

-x0,yl,b 

y:(rO),a 

do 

n2,_end_grp 


do 

r2,_end_bfy 


subl 

b,a x:(r0),b 

b,y:(r4) ;Radix 2 DIT butterfly 

kernel 



mac 

-xl,x0,b x:(r0)+,a 

a,y:(r5) 

macr 

-yl,y0,b x:(rl),xl 


subl 

b,a b,x:(r4)-»^ 

y:(r0),b 

mac 

xl,y0,b 

y:(rl)+,yl 

macr 

-x0,yl,b a,x:(r5)+ 

y:(rO),a 

_end_bfy 



move 

(rl)+nl 


subl 

b,a x:(r0',b 

b,y;(r4) 

mac 

-xl,xO,b x:(r0j>n0 

,a a,y;(r5) 

macr 

-yl.yO.b x;{rl),xl 

y:(r6),y0 

subl 

b,a b,x:(r4)+n4 y:(r0),b 

mac 

xl,y0,b xr(r6)+n6,x0 y:{rl)+,yl 

macr 

-x0,yl,b a,x:(r5)4n5 y:(r0),a 

_end_grp 



move 

nl,bl 
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Isr b n2,al ;divide butterflies per group by two 

Isl a bl,nl ;multiply groups per pass by two 

move al,n2 
_endjpass 

t 

; Do last IFFT pass 




move //2,n0 
move nO,nl 
move f/points/4,n4 
move n4,n5 
move «!>data,rO 
move tfodata,r4 
move r4,r2 
lua (rO)+,rl 
lua (r2)+n2,r5 
move C/0,m4 
move m4,m5 
move Ifcoef,r6 
move (r5)-n5 
move X:( 

move X:( 


;initialize pointer offsets 


rl) 

r5) 


output pointer A offset 
output pointer B offset 
initialize A inpi:t pointer 
initialize A output pointer 
save A output pointer 
initialize B input pointer 
initialize B output pointer 
bit-reversed addressing for 
bit-reversed addressing for 
initialize C input pointer 
predecrement output pointer 


,xl 

,a 


y:(r6),y0 

y:(rO),b 


output ptr. A 
output ptr. B 


do n2,_lastpass 

mac xl,yO,b x:{r6)+n6,x0 y:(rl)+nl,yl 
kernel 

macr -xO,yi,b a,x:(r5)+n5 y;(rO),a 

group 


;Radix 2 DIT butterfly 
;with one butterfly per 


subl b,a 
neg b 

move b,y:(r4) 
move x:(rO),b 
neg a 

move a,y;(r5) 
mac -xl,xO,b 
macr -yl,yO,b 
subl b,a 
lastpass 
move 
endm 


x:(rO)+nO,a 
x:(rl),xl 
b,x:(r4)+n4 

a,x:(r5)+n5 


y:(r6),y0 

y:(rO),b 


;complete last butterfly 
;and output conjugate 
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APPENDIX C. 


SINE COSINE GENERATOR 


This progr&m originally available on the Motorola DSP bulletin board. 
It is provided under a DISCLAIMER OF VARRANTY available from 
Motorola DSP Operation, 6501 Wm. Cannon Drive W., Austin, Tx., 78735. 

Sine-Cosine Table Generator for FFTs. 


Last Update 25 Nov 86 Version 1.2 


w 

sincos 

macro 

points,coef 

sincos 

ident 

1,2 



sincos 


macro to generate sine and cosine coefficient 
lookup tables for Decimation in Time FFT 
twiddle factors. 


points 


number of points (2 - 32768, power of 2) 


coef 


base address of sine/cosine table 
negative cosine value in X memory 
negative sine value in Y memory 


Latest revision - 25-Nov-86 


pi 

equ 

3.141592654 

freq 

equ 

2.0*pi/@cvf(points) 


org 

x:coef 

count 

set 

0 


dup 

points/2 


dc 

-0COS(@cvf(count)* freq) 

count 

set 

endm 

count+1 


org 

y:coef 

count 

set 

0 


dup 

points/2 


dc 

-esin(0cvf(count)*freq) 

count 

set 

endm 

counts1 


endm 

;end of sincos macro 
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APPENDIX D. 


16 QAH DATA FILES 


The following two data files are prescaled (full scale constellation 
values divided by 256) and conjugated. Both data files must be loaded 
into the DSP prior to running QAM16ENC. They are properly formatted to be 
read into memory. They occupy the first 16 locations of x and y memory. 


Data file '16QAMRE.DAT' 

HX 

00000000,OOOOOOOF 
FFF800,FFE800,FFF800,FFE800 
000800,000800,001800,001800 
FFF800,FFF800,FFE800,FFE800 
000800,001800,000800,001800 


Data file '16QAMIM.DAT' 

HY 

00000000,OOOOOOOF 
000800,000800,001800,001800 
000800,001800,000800,001800 
FFF800,FFE800,FFF800,FFE800 
FFF800,FFF800,FFE800.FFE800 
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APPENDIX E. 


16 QAM ENCODER 


This program is an encoder for a 16 QAM system. A symbol consisting of 
four bits is read into the engwder and using a look up table outputs the 
values for the points on the constellation and are loaded into an array. 
The data files 16QAMRE.DAT and 16QAM1M.DAT must be loaded prior to 
running 16QAMEN. Once the 256 values are loaded into an array an IFFT 
is performed to generate the MFH signal. 

The following registers are modified: 
r2 
r3 
r4 
r5 

The following Data ALU Registers are modified: 


aO 

al 

a 

xO 

xl 

yo 

bO 

bi 

b 


qaml6en 

ident 1,1 
page 132,5A 

opt nomd,nomex,loc 

include 'sincos' ;include macro for sine cosine values 
include ’ifft' ;include macro for calculating IFFT 


; Define memory locations to be used in the program 


start 

equ 


$100 

;starting program location 

startifft 

equ 

$350 

;starting program location for IFFT 

input 

equ 


40 

;starting location of input bits 

output 

equ 


50 

;starting location of output bits 

coef 

equ 


128 

;location of coefficients for IFFT 

points 

equ 


256 

;number of points for IFFT 

data 

equ 


768 

;location of input data for IFFT 

odata 

equ 


1280 

;location of output data for IFFT 

locate 

equ 


$9C4 

;location of message to be sent (hex) 

sincos 

points,coef 

opt 


mex 
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org pjstart 

move //locate, r3 

begin 

move //data,r4 

do #201,cod 

move //output, r5 

jsr readbit 

move r4,x0 

move H/data+lDOja 

cmp x0,a ;check to see if 100 symbols have been read 

jseq outzero ;output 55 zeros, read remaining symbols 

jsr outbit 

code 

org p:startifft 

ifft points,data,odata,coef ;macro call to perform IFFT 
; jmp begin 

nop 

swi 


; The subroutine readbit reads the input value into memory one bit at a 
; time 


readbit 


loop 


move 

//input+3,r2 

move 


move 

//>$l,x0 

do 

//4, loop 

and 

x0,a a 

move 

al,x:(r2)- 

move 

xl ,a 

asr 

a 

rts 



y:(r3)+,a ;read input symbol 
; symbol is read in one bit at a time 


; The subroutine outbit outputs the values stored by readbit. This four 
; bit symbol is used as an index into a lookup table. 


outbit 

move 

//input, r2 


clr 

b 


clr 

a 


addl 

b,a 


do 

//4, loop 2 

;symbol 

move 

x:(r2)+,b0 


addl 

b,a 


loop2 

move 

a0,r5 

;symbol 

nop 

move 

x: (r5),x0 


move 


y:(r5),y0 


read out one bit at a time 

is used to index into lookup table 

;get real constellation value 
:get imaginary constellation value 


39 




move 

move 

rts 

xO,x:(r4) 

;load real value into IFFT array 
y0,y:(r4)4 ;load Imag value into IFFT array np 

outzero 

do 

//SS, endzero 

;subroutine loads zeros in center 

move 

»0,x0 

;55 locations of IFFT array to 

move 

nop 

move 

endzero 

nop 

rts 

xO,x:(r4) 

;allow for filtering 

x0,y;{r4)+ 
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APPENDIX F. 16 QAM BOUNDARY FILE 


The following data file BOUN16.d contains the sjnnbol data for each bounded 
area used in the 16 QAM decoder. It must be loaded into the DSP prior to 
running 16QAM0EC. 

Data file 'B0UN16.D’ 

HX 

00000200,0000020? 

000000,000008,000000,000004 
OOOOOD,OOOOOA,000001,000006 
OOOOOE,000009,000002,000005 
00000F,00000B,000003,000007 
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APPENDIX G. 16 QAM DECODER 


This program is a decoder for 16 QAM system. After computing the FFT 
of the received data the points are read out of the array one at a time 
and compared with boundaries on the constellation to determine the 
received point. The decoded symbols are then stored in memory. It is 
assumed that the sampled input data is loaded into the input data 
location (1280) of the FFT. If the sampled data is stored elsewhere a 
routine to load the FFT array will need to be added. 

The following registers are modified: 

rO r4 

rl r5 

r2 r6 

The following Data ALU registers are modified: 

a b 

xO xl 

yO yl 


qaml6dec 

ident 1,1 

page 132,66,3,3,0 

opt nomd,nomex,loc,nocex,mu,cex 


include 

include 


'sincos' 

'fff 


org 1:$0000 

location dsm 16 
input dsm 16 

endlong equ ' 
org x:endlong 
storr6 ds 1 
org x:512 


;include macro for sine cosine values 
;include macro for calculating FFT 


;storage locations for look up table 
;storage locations for input data 


boundaryl 

ds 

4 

;reserves 

four 

locations 

for boundary pts. 

boundary2 

ds 

4 

;reserves 

four 

locations 

for boundary pts. 

boundary3 

ds 

4 

;reserves 

four 

locations 

for boundary pts. 

boundary^ 

ds 

4 

;reserves 

four 

locations 

for boundary pts. 


startfft equ $100 


:starting location for fft program 




start equ $250 ;starting location for decoder program 
points equ 256 ;number of points for FFT 

coef equ 1024 ;location of coefficients for FFT 

data equ 1280 ;location of input data for FFT 

odata equ 1536 {location of output data for FFT 

{Full scale values for 16 QAM constellation. 

four equ $200000 

three equ $180000 

two equ $100000 

one equ $080000 

zero equ $000000 

mone equ $F80000 

mtvo equ $F00000 

mthree equ $e80000 

mfour equ $e00000 

sincos points,coef 
opt mex 

org p:startfft 

fft points,data,odata,coef {Macro call for FFT 

opt cex 
org p:start 
jsr initialize 
do //201,_endrun 
jsr readdata 
jsr findbound 
_endrun 

nop 

swi 

{This initialization routine initializes register and modifiers 
{as well as clearing the memeory. 

{The constellation is also loaded into memory here. 

initialize 

move #$ffff,m0 {reset register to linear addressing 

move <!^$ffff,ml {reset register to linear addressing 

move i^$ffff,m2 {reset register to linear addressing 

move #$ffff,m4 {reset register to linear addressing 

move (/$ffff,m5 {reset register to linear addressing 

move #15,m6 {set register for modulo 15 addressing 

move #0,rl 

clr b #$0,r0 

clr a r0,r5 

do #50,clrmem 

move a,x:(r0)+ b,y:(r5)+ {dear first 50 memory locations 
clrmem 
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DOve f/inputibl 

move bl,x:storr6 

Now load full scale values of the constellation in the table locations. 


move 

/^location, 

rO ;Real 

Imag 

move 

r0,r4 




move 

#mone,a 




move 

tfmone,b 




move 

a,x:(rO)+ 

b,y:(r4)+ 

;-l 

-1 

move 

(!tmthree,b 




move 

b,x: (roi-f 

a,y:(r4)+ 

;-3 

-1 

move 

a,x:(rO)+ 

b,y:(r4)+ 

;-l 

-3 

move 

#mthree,a 




move 

a,x:(rO)+ 

b,y;(r4)+ 

;-3 

-3 

move 

{t>mone,b 




move 

Itone, a 




move 

a,x:(r0)+ 

b,y:(r4)-f 

; 1 

-1 

move 

li^mthree.b 




move 

a,x:(r0)+ 

b,y:(r4)+ 

; 1 

-3 

move 

f!^three,b 




move 

//mone,a 




move 

b,x:(r0)+ 

a,y:{r4)+ 

; 3 

-1 

move 

//mthree,a 




move 

b,x:(r0)+ 

a,y:(r4)4 

; 3 

-3 

move 

#mthree,a 




move 

#three,b 




move 

#mone,xO 




move 

ffone.yl 




move 

xO,x:(r0)+ 

yl,y:(r4)+ 

-1 

1 

move 

xO,x;(r0)+ 

b,y;(r4)+ 

-1 

3 

move 

a,x:{r0)+ 

y/,y:(r4)+ 

-3 

1 

move 

a,x:(r0)4 

b,y:(r4)+ 

-3 

3 

move 

//one, a 




move 

//one,xO 




move 

//three, yl 




move 

xO,x:(r0)+ 

a,y:(r4)+ 

1 

1 

move 

b,x:(r0)4 

a,y:(r4)+ 

3 

1 

move 

xO,x: (rO)-t- 

yl,y:(r4)+ 

1 

3 

move 

b,x:(r0)4 

yi.yj(r4)+ 

3 

3 

move 

jj/odata,xO 




move 

tf$eff,yO 




move 

x0,y:$424 

;store location 

for input to decoder 

move 

y0,y:$425 

;store location 

for output of decoder 

rts 






The subroutine readdata reads in the data from the output of the FFT. 
The values are read out one point at a time. The values are coiq>ared to 
boundaries on a partioned constellation to determine the received point. 
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readdata 


move 

y:$424,r0 


nop 



move 

rO,yO 


move 

//odata-flOOid 

;check to see if first 100 values have been 

read 



cmp 

yO,a 


jseq 

_delzero 

;delete 55 zeros 

move 

x:storr6,r6 


move 

x:(rO),a 

;input data is loaded into storage location 

move 

y:(rO)+,b 

:for decoding 

move 

a,x:(r6) 


move 

b,y: (r6) + 


move 

r6,x:storr6 


move 

r0,y:$424 


rts 

• 



> 

;subroutine 

• 

_delzero removes the zeros installed in the encoder 

9 

_delzero 



move 

#odata+155,r0 


nop 



rts 



9 

; The subroutine findbound 

compares the value that has been read out of 

; the array 

to the boundries on the constellation to decode the point. 

; First the magnitude alone 

is used to find the correct bounded area then 


the signed values are used to determine the correct quadrant is used to 
increment the boundary pointer to find the correct point. 


findbound 

move 

x:-{r6),a 

; real value is stored in a 

move 

//boundaryl,r2 ; load starting position for bounded values 

move 

tftwo.yO 


cmpm 

y0,a 

y:(r6),b ; compare mag of real value to two 

jgt 

bigtvo 

; imaginary value is stored in b 
; x>2 

cmpm 

yO.b 

x:(r2),x0 ; compare mag of imag value to two 

jit 

continue 

; x<2,y<2, load r2 with boundary 1 and 

move 

(!>four,xl 

; continue 

cmpm 

xl,b 

f/boundary3,r2 ; compare magnitude of imag value 

nop 

move 

x: (r2),x0 


jmp 

continue 

;x<2,y>2 and y<4, load r2 with boundary4 

bigtvo 

cmpm 

y0,b 

//boundary2,r2 ;x>2 y<2 


nop 
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move x:(r2),x0 

jit continue ;x>2 and x<4, y<2 load boundary2 and 

continue 

cmpm xl,b #boundary4,r2 

nop 

move x:(r2),x0 

jmp continue ;x>2»y<4 load boundary 4 and continue 

; This part of the routine finds the correct quadrant and updates the 
; pointer to the correct point. 

continue 

clr a x:(r6),xl 

cmp xl,a y:(r6),yl 

jgt negx 

cmp yl,a //3,n2 

jgt posxnegy 
posxposy 

jmp outputdata ;output is in first quadrant 

posxnegy 

move x: (r2)-t-n2,x0 ;update r2 by 3, fourth quadrant 
jmp outputdata 

negx 

cmp yl,a iH,n2 

jgt negxnegy 
negxposy 

move x:(r2)-fn2,x0 ;update r2 by 1, second quadrant 
jmp outputdata 
negxnegy 

move //2,n2 
nop 

move x:(r2)+n2,x0 ;update r2 by 2, third quadrant 

outputdata 

move y:$425,r0 

move x:(r2),a 

nop 

move al,y:(rO)> :move decoded symbol to output location 

move r0,y:$425 

rts 
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APPENDIX H. 


2/3 RATE CODE DATA FILES 


The following two data files are prescaled (full scale constellation 
values divided by 256) and conjugated. Both data files must be loaded 
into the DSP prior to running 23ENC0D. They are in the proper format to 
be loaded into memory. 


Data file '23REAL.DAT' 

HX 

00000000,00000007 
001000,FFFOOO,000800,FFF800 
FFF800,000800,FFFOOO,001000 


Data file ’23IMAG.DAT’ 

HY 

00000000,00000007 

FFF800,000800,FFFOOO,001000 

FFFOOO,001000,FFF800,000800 
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APPENDIX I. 2/3 RATE CODE ENCODER 


This program is an encoder for a 2/3 rate convolutional code. A symbol 
consisting of two bits is read into the encoder, combinational 
logic is used to generate the third bit then a look up table is used 
to output the values from a point on the constellation. The data files 
23REAL.DAT and 23IMAG.DAT must first be loaded in memory prior to 
23ENC0D. The 256 values are loaded into an array and an IFFT is 
performed to generate the MFM signal. 

The following registers are modified; 
rO r3 

rl r4 

r2 r5 

The following Data ALU registers are modified; 
a aO al 
b bO 
xO xl 
yO yl 


23encod 

ident 1,1 
page 132,54 

opt nomd,nomex,loc 


include ’sincos' 
include 'ifft' 


start 

equ 

$100 

;starting program location 

startifft 

equ $350 

;starting program location for IFFT 

points 

equ 

256 

;number of points for IFFT 

coef 

equ 

128 

;location of coefficients for IFFT 

data 

equ 

768 

;location of input data for IFFT 

odata 

equ 

1280 

;location of output data for IFFT 

locate 

equ 

$9C4 

;starting location for message to be sent 

input 

equ 

40 

;storage locations for input bits 

output 

equ 

50 

;storage locations for encoded bits 

statemeu: 

equ 60 

;storage for past bits 


sincos points,coef 
opt mex 


org p:start 

move ^/locate, r3 
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begin 

move //data,r4 
do f/201,code 

move //output, r5 

jar readbit 

jsr encode 

move r4,x0 

move //data‘4-100,a 

cmp xO,a ;check to see if 100 symbols have been read 

Jseq outzero ;output 55 zeros, read remaining symbols 

jsr outbit 

code 

org p:startifft 

ifft points,data,odata,coef;macro call to perform IFFT 

; jmp begin 


nop 

svi 

« 

; The subroutine readbit reads two bits from the message and stores them 

; in memory. 

• 

9 

readbit 
move 
move 
move 
do 
and 
move 
move 
asr 

loop 

rts 

9 

; The subroutine encode performs the convolutional encoding to generate 

; a redundant bit. 

• 

9 

encode 


move 

//input, rO 


move 

//output, r5 


move 

//statemem,rl 


move 

x: (rO)-f ,x0 

;read bit xO 

move 

x:(r0)-,xl 

;read bit xl 

move 

x:(rl)4,a 

;read past bit si 

move 

x:(rl)-,b 

;read past bit s2 

eor 

xl,a 

;xl eor si 

move 

a,yl 


eor 

yl,b 

;xl eor si eor s2 = y2 

move 

x0,y:{r5) 

;store xO 

move 

x;(rl)+,a 

;read past bit si 

move 

x:(rl)-,x0 

;read past bit s2 


i//input+l,r2 

y;(r3)+,a 

//>$l,x0 
//2, loop 

xC,a a,xl 

al,x:(r2)> 
xl,a 
a 
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eor 

xl,a 

;xl eor si » yl 

move 

x0,x:(rl) + 

;move s2 to si 

move 

xl,x:(rl) 

;move xl to si 

move 

a,x: (rO)-f 

:temp store yl 

move 

b,x:(rO) 

;temp store y2 

rts 




; The subroutine outbit reads the three bits out of memory to form a 
; symbol. The symbols are reconstructed by reading the bits one at a 
; time added to to an empty register and shifted left until three bit 
; symbol is formedA lookup table is used to get the values for the point 
: cn the constellation. 

outbit move #input,r2 

clr b 

clr a y:(r5),b0 

addl b,a 

do //2,loop2 

move x:(r2)+,b0 

addl b,a 

loop 2 

move a0,r5 ;move symbol to r5 

nop 

;use value in r5 to get points off the constellation 

move x:{r5),x0 

move yj(r5),y0 

;move constellation points to IFFT array 

move x0,x:(r4) 

move y0,y:{r4)+ 

nop 
rts 
J 

;subroutine outzero puts 55 zeros in IFFT array for filtering 


outzero 

do //55,endzero 

move <lfO,xO 

move x0,x:(r4) 

nop 

move x0,y:{r4)+ 

endzero 
nop 
rts 
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APPENDIX J. 2/3 RATE CODE BOUNDARY FILE 


The following data file B0UN23.0 contains the points on the constellation 
used to find the minimum distances to states after the proper quadrant is 
found. 

Data file 'BOUN23.D 
HX 

00000200,0000020F 

000040,0000A2,000044,000047 

000041,000042,000044,000046 

000041,000043,000045,000046 

000040,000043,000045,000047 
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APPENDIX K. 2/3 RATE CODE VITERBl DECODER 


; This program is a Viterbi Decoder for a 2/3 rate convolutional encoder. 
; There is a 16 time period delay which will approach the maximum possible 
: gain for this type of encoder. There are 64 locations needed in memory 

; (16 past time periods x 4 states s 64). 

• 

> 

; The following registers are modified: 

: rO rS 

; rl r6 

; r2 r7 

; t4 

• 

f 

; The following Data ALU registers are modified: 

; a b xO yO 

; al bl xl yl 

• 

I 

; The following register modifiers are used; 

; ml nO 

; m5 n2 

; m6 

• 

f 

mfm23dec 

ident 1,1 

page 132,66,3,3,0 

opt nomd,nomex,loc,nocex,mu,cex 

include ’sincos' 
include ’fft' 


org 1:$0000 


period dsm 

64 

:64 storage locations 

location dsm 

8 

;constellation points 

input dsm 

16 

;past 16 input bits 

tables dsm 

4 

;accumlated distance 

temp dsm 

4 

;temp storage for distance table 

endlong equ 

* 


org x:endlong 

storr6 ds 

1 


ynow dsm 

3 

;input bits 

org y:endlong 

ypast dsm 

2 

;past bits 

org x:512 

boundary! ds 

16 

;storage for boundary data 

startfft equ 

$100 

;starting location for FFT 
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points 

equ 

256 

;number of points for FFT 


coef 

equ 

1024 

;location of coefficients for 

FFT 

data 

equ 

1280 

:location of sampled data for 

FFT 

odata 

equ 

1536 

;output location for FFT 


start 

equ 

$250 

{Starting location for decoder 

program 

; Define 

full scale constellation values. 


two 

equ 

$100000 



one 

equ 

$080000 



zero 

equ 

$000000 



mone 

equ 

$F80000 



mtvo 

equ 

$FOOOOO 



large 

equ 

.9 



small 

equ 

.1 



offset 

equ 

$000000 

;can be used to distort data 



sincos points,coef 
opt mex 

org prstartfft 
fft points,data,odata,coef 
opt cex 
org p:start 
jsr initialize 
do #201,_endrun 
jsr readdata 
jsr findmindist 
jsr accumdist 
jsr traceback 
jsr outputdata 
_endrun 

nop 

svi 

f 

;this initialization routine initializes register and modifiers 
;as well as clearing the memeory. The constellation is also loaded 
;into memory here. The accumulated distance array is set so that 
;state zero starts out at a value of zero and all others start out 
;larger, forcing the paths to merge at the zero states. 

> 

initialize 


move 

j|^$ffff ,m0 

;sets linear addressing 

move 

P63,ml 

;sets modulo 63 addressing 

move 

#$ffff,m2 

;sets linear addressing 

move 

#$ffff,m4 

;sets linear addressing 

move 

//63,m5 

;sets modulo 63 addressing 

move 

n5,m6 

;sets modulo 16 addressing 

move 

//0,rl 
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i 


clr b #$0,r0 

cli* a rO,r5 

do 1/256,clrmem 
move a,x:(rO)+ b,y:(r5)+ 

clrmem 

move #tables-»'l,r7 

move #$400000,al 

rep #3 

move al,x:(r7)+ 

move #input,bl 

move bl,x:storr6 

move #odata'*16,r0 

move r0,r2 

move #0,x0 

move #0,y0 

do #16,_clrreg 

move x0,x:{r0)+ 

move y0,y:(r2)+ 

_clrreg 

;Now load full scale values of the constellation in the table locations. 


move 

#location. 

rO ;Real 

Imag 

move 

r0,r4 




move 

#two,a 




move 

if one , yl 




move 

#mtwo,b 




move 

#mone,yO 




move 

#one,xl 




move 

#mone,xO 




move 

a,x;(r0)+ 

yl,y:(r4)+ 

2 

1 

move 

xl,x:{r0)+ 

a,y:(r4)+ 

-2 

-1 

move 

xO,x:(r0)+ 

a,y:(r4)+ 

1 

2 

move 

b,x:(r0)+ 

yl,y:(r4)+ 

-1 

-2 

move 

b,x:(r0)+ 

y0,y:(r4)+ 

;-i 

2 

move 

xl,x:(r0)'f 

b,y:(r4)+ 

1 

-2 

move 

x0,x:(r0)4 

b,y:(r4)+ 

1 

-2 

move 

a,x: (r0)-»^ 

y0,y:(r4)+ 

2 

-1 

move 

#odata,x0 




move 

#$eff,y0 




move 

x0,y:$424 




move 

y0,y:$425 





rts 

readdata reads in the data from the outpu of the FFT. The data is read 
in as complex points on the constellation. 
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readdata 


move 

nop 

y:$424,r0 

move 

r0,y0 

move 

//odata-fl00,a 

cmp 

y0,a 

jseq 

_delzero 

move 

x:storr6,r6 

move 

(I^>offset,x0 

move 

X:(rO), a 

add 

x0,a y 

add 

x0,b a 

move 

b,yj(r6)+ 

move 

r6,x:storr6 

move 

rts 

r0,y:$424 


.delzero 

move 

nop 

rts 


#odata4l55,r0 


the minimum distance is found to the closest point in every state and 
stored. The values are stored so that indexing is made easier, state 
0,2,3,1. This will greatly reduce the number of cycles needed later, 
a smoothing function is used to accumulate distances in the accumulated 
table so this minimum distance is multiplied by .1. 

The subroutine findmindist finds the quadrant of the received point 
which is used as a pointer into the boundary table to read the four 
closest points, one in each state. 

findmindist 


;x<0 


move 

#boundaryl,r2 

clr 

a x:-(r6),xl 

cmp 

xl,a y:(r6) 

jgt 

negx 

cmp 

yl,a //12,n2 

jgt 

posxnegy 

posxposy 

jmp 

findist 

posxnegy 

move 

x:(r2)-fn2,x0 

jmp 

findist 

negx 

cmp 

yl,a #4,n2 

jgt 

negxnegy 

negxposy 

move 

x: (r2)-t-n2,x0 

jmp 

findist 
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A 


negxnegy 

move //8,n2 
nop 

move x:(r2)+n2,x0 

« 

> 

; findist finds the distance from the received point to the points read 

; from the boundary table. 

• 

findist 


move 

x:(r2)-f ,r0 

move 

#tables,r4 

move 

x:(rO),a 


sub 

xl,a 

y:{rO),b 

sub 

yi.b 

a,xO 

mpy 

xO,xO,a 

b,yO 

mac 

yO.yO.a 

x: (r2) + ,r0 

move 

//small, xO 

a.yO 

mpy 

xO,yO,a 


move 

x:(rO),a 

a,y:(r4)+ 

sub 

xl ,a 

y:(rO),b 

sub 

yi.b 

a,xO y:(r 

mpy 

xO,xO,a 

b,yO 

mac 

yO.yO.a 

x: (r2) + ,r0 

move 

//small,xO 

a.yO 

mpy 

xO,yO,a 


move 

x:(rO),a 

a,y;(r4)+ 

sub 

xl,a 

y:(rO),b 

sub 

yi.b 

a,xO 

mpy 

xO,xO,a 

b,yO 

mac 

yO.yO.a 

x: (r2) + ,r0 

move 

//small, xO 

a.yO 

mpy 

xO,yO,a 


move 

x:(rO),a 

a,y:(r4)- 

sub 

xl,a 

y:(rO),b 

sub 

yi.b 

a,xOy;(r4)-,y0 

mpy 

xO,xO,a 

b.yO 

mac 

yO,yO,a 

x: (r2) + ,r0 

move 

//small,xO 

a.yO 

mpy 

xO,yO,a 


move 

rts 


a,y;(r4) 


;the accumulted distance routine adds the smallest distance from the 
;prevlously computed table for all pathes going into a state and 
;does this for all four states. Since only certain transitions are 
;alloved the calculations are done in a specific order to reduce delay. 

f 

accumdist 

clr a /^tables.rO 
move l!f$7fffff,al 

move r0,r4 
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move 

tftemp,r2 

move 

in,mO 

move 

m0,m4 

move 

in,nl 

move 

nl,n5 

move 

rl,r5 


;find minimum distance 

to state zero 

do 

)!t2, statezero 


move 

x;(rO),xO 

y:(r4),b 

add 

xO,b 


cmp 

b,a 


tge 

b,a 

r0,r3 

tge 

b,a 

r4,r7 

move 

x:{rO)+,xO 

y:(r4)+,b 

statezero 



move 

r3,x:(rl)+nl 


move 

a,x: (r2)-f 

y:(r4)+,b 

clr 

a 

r7,y:(r5)+n5 

move 

#$7fffff,al 



;find minimum distance 

to state two 

do 

//2, statetwo 


move 

x:(rO),x0 

y:(r4),b 

add 

xO,b 


cmp 

b,a 


tge 

b,a 

r0,r3 

tge 

b,a 

r4,r7 

move 

x: (rO) + ,xO 

y:(r4)+,b 

statetwo 



move 

r3,x:(rl)+nl 


move 

a,x: (r2)-t- 

y: (r4) + ,b 

clr 

a 

r7,y:(r5)+n5 

move 


#tables+2,r4 

move 


r4,r0 

move 


x: (rl)--nl,a 

clr 

a 

x:(rl)-,b 

move 

/;$7fffff,al 


move 

rl,r5 


;find minimum distance 

to state one 

do 

(!^2,stateone 


move 

x; (rO),x0 

y:(r4),b 

add 

xO,b 


cmp 

b,a 


tge 

b,a r0,r3 

tge 

b,a r4,r7 
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move x: (rO) + ,xO y;(r4)-,b 
stateone 

move r3,x:(rl)+nl 
move a,x:(r2)-»- 

clr a r7,y:(r5)+n5 

move tf$7fffff,al 

;find minimum distance to state three 
do #2,statethree 
move x:(rO),x0 y:(i4),b 

add xO,b 
cmp b,a 

tge b,8 r0,r3 

tge b,a r4,r7 

move x:(rO)+,xQ y:{r4)-,b 

statethree 

move r3,x:{rl)+nl 
move a,x:(r2)4 

clr a r7,y:(r5)+n5 

move j'/$7fffff,al 
move (r4)+ 

;now move new accumulated distances into the accumulated distance 
;table from the temporary table 

;also find the min distance state and store in r4 which is no longer used 


move 


#$ffff,mO 

move 


#$ffff,m4 

move 


if temp , r3 

move 


stables,rO 

move 


//large, xl 

move 


#2,n0 

do 


ll/3,endtable 

move 


x: (r3) + ,x0 

mpy 


xl,xO,a 

cmp 

a,b 

a,x:(rO)+nO 

tge 

a,b 

r0,r4 

endtable 

move 


//tables-t-l,rO 

do 


#4,endtablex 

move 


x;(r3)+,x0 

mpy 


xl,xO,a 

cmp 


a,b a,x:(r0)4n0 

tge 


a,b rC,r4 

endtablex 

;store in 

rO 

instead of r4 


move r4,r0 

move tf4,nl 
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move 

rts 


(rO)-nO 


;the traceback routine now goes back through every time period starting 
;vith the current time period and finds the state from which the path 
;came from one time period previous. At the end of this search, the 
;last state found will also point to the path at that state, which is the 
;output of the trellis. 

traceback 

;find the displacement from the pointer to table and store value in n4 

move /i^tables,nO 

move (rl)-nl 

lua (r0)-n0,n5 

move rl,r5 

do //15,endtrace 

move (rl)-nl 

move x:(r5+n5),r0 

move rl,r5 

lua (rO)-nO,n5 

endtrace 

move //location, rO 

move y:(r5+n5),a 

rts 

;the output data routine unscrambles the path order and finds one 

;of the two points on the constellation coresponding to the output state 

;which is closest to the original input at that time period. 

outputdata 

move a,b 

move //>$bl,xO 
cmp xO,a //>$b2,y0 

teq yO,b 

cmp yO,a //>$b3,x0 

teq xO,b 

cmp xO,a //>$bl,yO 

teq yO,b 

move //>$b5,x0 
cmp xO,a f/>$b7,y0 

teq yO,b 

cmp yO,a 

teq xO,b 

move b,r2 
move //tables,n2 
move x:storr6,r6 
lua (r2)>n2,n3 
move n3,a 
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asl 

a 

asl 

a 

move 

a,n0 

move 

r6,r3 

lua 

(r0)+n0,r4 

move 

#>$7fffff,xl 

move 

r4,r0 

do 

//4, endout 

move 

.'';{r3),a y:(r6) 

move 

x:(r0)+.x0 y:(r4) 

sub 

x0,a 

sub 

y0,b a,x0 

mpy 

x0,x0,a b,y0 

mac 

y0,y0,a 

tfr 

a,b xl,a 

cmp 

xl,b 

tit 

b^a r0,r7 

move 

a,xl 

endout 


clr 

a (r7)- 

move 

#location,n0 

move 

r7,r0 

move 

•#$f,al 

lua 

(r0)-n0,r7 

move 

r7,x0 

and 

x0,a 

jsr 

convoldec 

move 

y:$425,r0 

nop 


move 

a0,y:(r05+ 

move 

r0,y:$425 

rts 



;The subroutine convoldec decodes the received symbol by using 
;combinational logic. 


convoldec 


move 

# 3 mow+ 2 ,r 0 

move 


move 

#ypast,r7 

do 

#3,loop 

and 

x0,a 

move 

al,x:(rO)- 

move 

xl,a 

asr 

a 

loop 



move x:(rO)+,a y:{r7)+,y0 
move x:(rO)-f,b y:(r7)“,yl 
move x:(rO)-,a 


;read si 
;read yO and s2 
;read yl 
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nop 




eor 

yO,a 


;yl eor si =xl 

move 


yl,y:(r7)+ 

;update past states 

move 

a,y:(r7)- 



move 

b,x:(r0)“ 



move 

a,x:(rOi+ 



clr 

a 



clr 

b 



do 

#2,loop2 



move 

x: (rO)-,bO 


;catput decoded bits 

addl 

b,a 




loop2 

rts 









APPENDIX L. V.32 DATA FILES 


The following data files are prescaled (full scale constellation values 
divided by 256) and conjugated. Both data files must be loaded into the 
DSP prior to running MFMENC. They are properly formatted to be read into 
memory. 

Data file •QAMREAL.DAT' 

HX 

00000000,OOOOOOIF 
FFEOOO,000000,000000,002000 
002000,000000,000000,FFEOOO 
FFF000,FFF000,001000,001000 
001000,001000,FFF000,FFFOOO 
FFE800,000800,FFE800,000800 
001800,FFF800,001800,FFF800 
000800,FFE800,000800,000800 
FFF800,001800,FFF800,FFF800 

Data file ’QAMIMAG.DAT’ 

HY 

00000000,OOOOOOIF 
FFF800,001800,FFF800,FFF800 
000800,FFE800,000800,000800 
FFE800,000800,FFS800,000800 
001800,FFF800,001800,FFF800 
001000,001000,FFFOOO,FFFOOO 
FFF000,FFF000,001000,001000 
FFEOOO,000000,000000,002000 
002000,000000,OOCOOO,FFEOOO 
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APPENDIX H. 


V.32 ENCODER 


; This is a an encoder for V.32 standard with differential encoding. 

; The encoder can also be used with the differential removed by commenting 
; out the call to diff. 

mfmencod 

ident 1,1 


page 

132,54 

opt 

nomd,nomex,loc 

include 

'sincos* 

include 

'ifft' 


org x:$40 
statemem ds 

3 

;set 3 locations for p >.st states 

input 

ds 

4 

;set up 4 locations for input bits 

org 

ylpast 

y:$40 

ds 

2 

;set up two locations for diff encoder 

output 

ds 

1 

;output of convolutional encoder 

start 

equ 

$100 

;starting location for encoder 

startifft equ 

$350 

;starting location for IFFT 

points 

equ 

256 

;number of points for IFFT 

odata 

equ 

1280 

;output of IFFT 

data 

equ 

768 

;input data for IFFT 

coef 

equ 

128 

;location for coeficients for IFFT 

locate 

equ 

$9C4 

;starting location for message 


sincos points,coef 
opt mex 


org 

p:start 

move 

}/ylpast,r5 

move 

f/statemem,r3 

move 

{!^locate,r6 


begin 


move 

l!>data,r7 

do 

l!<201,code 

move 

l!>output,r4 

jsr 

readbit 

jsr 

diff 

jsr 

encode 


;reads 201 message symbols 
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move r7,x0 

move //data4l00,a 

cmp xO, a 

jseq outzero 

jar outbit 

code 

org p;startifft 

ifft points,data,odata,coef 

: jmp begin 

nop 

svi 


jchecks if first 100 symbols read 
{Outputs 55 zeros 


{subroutine readblt reads in the four bit symbol one bit at a time 
readbit 

move tfinput+3,r2 

move y:(r6)+,a 

move ('/>$!,xO 

do , loop 

and x0,a a,xl 

move al,x:(r2)- 

move xl,a 

asr a 

loop 

rts 

t 

{subroutine diff differentially encodes the two most significant bits 


diff 


move 

/7input,rl 


move 


y:{r5)+,y0 

move 

x;(rl)+,a 

y:(r5)-,yl 

move 

x;(rl)-,b 


eor 

y0,a 

a,x0 

eor 

yi,b 

a,xl 

move 


x0,a 

and 

y0,a 

b,yl 

eor 

yl,a 

xl,b 

move 

b,x: (rl)-f 

b,y:(r5)+ 

move 

a,x:(rl) 

a,y;(r5)- 

rts 




subroutine encode convolutionally encodes fhe four bits to generate a 
fifth bit 


encode 

move »'/input,rO 
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w 

« # 


move 

output,r4 


move 

//statemen, 

rl 

move 

x: (rO)'f ,xl 


move 

move 

x:(rl )4,a 

a,y:(r4) 

and 

xl,a 

x: (rO),x0 

move 

x: (rl)-,b 


eor 

xO,b 

a,yO 

eor 

yO.b 

b,yl 

move 

b,x:(rl)4 

y:(r4),b 

and 

yi.b 

xO,a 

move 

(rl)4 


eor 

xl,a 

x:(rl),x0 

eor 

xO,a 

y:(r4),yl 

move 

b,yO 


eor 

yO,a 

yltx:(rl)- 

move 

rts 

a,x:(rl)+ 



subroutine outbit reads the five bit symbol and uses it to index into the 
lookup table to get the values of the point on the constellation 


outbit 

move //input, r2 

clr 

b 

clr 

a y:(r4),l 

addl 

b,a 

do 

//4, loop2 

move 

x;(r2)+,b0 

addl 

b,a 

loop2 

move 

a0,r4 

nop 

move 

x:(r4),x0 

move 

y: (r4),y0 


move the values read from the constellation to the input array for the 
IFFT 

move x0,x:(r7) 

move y0,y:{r7)+ 

nop 
rts 

subroutine outzero loads 55 zeros into the IFFT array for filtering 


outzero 

do f!>55,endzero 
move (^0,x0 
move x0,x:(r7) 
nop 
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move 

endzero 

nop 

rts 


xO, 


:{r7)+ 
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APPENDIX N. V.32 BOUNDARY FILE 


This data file BOUND.D contains the eight closest points (one in each 
state) to each of the 52 bounded areas used in the partitioned 
constellation. It must be loaded prior to using MFNOIFDE. 

Data file 'BOUND.D' 

HX 

00000200,0000039F 

000082,000086,00008B,00008D’ 

000093,C00095,00009A,00009E 
000082,000086,000089,00008F 
000093,000095,00009A,00009E 
000082,000086,000089,C0008F 
000091,000097,00009A,00009E 
000082,000086,00008B,00008D 
000091,000097,00009A,00009E 
000082,000086,00008B,00008D 
000093,000094,00009A,00009D 
000082,000086,000089,00008F 
000092,000095,000099,00009E 
000082,000086,000089,00008F 
000090,000097,000099,00009E 
000082,000086,00008B,00008D 
000091,000096,00009a,000090 
000083,000084,000088,000080 
000093,000094,00009a,000090 
000080,000087,000089,000Q8F 
000092,000095,000099,00009E 
000080,000087,000089,00008F 
000090,000097,000099,00009E 
000083,000084,00008B,00008D 
000091,000096,00009A,00009D 
000082,000085,00008A,00008D 
000093,000095,00009A,00009E 
000082,000085,000088,00008F 
000093,000095,00009A,00009E 
000081,000086,000089,00008E 
000091,000097,00009A,00009E 
000081,000086,00008B,00008C 
000091,000097,00009A,00009E 
000082,000085,00008A,00008D 
000093,000094,00009A,00009D 
000082,000085,000088,00008F 
000092,000095,000099,00009E 
000081,000086,000089,00008E 
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000090,000097,000099,00009E 
000081,000086,00008B,00008C 
000091,000096,00009A,00009D 
000082,000085,00008a,000080 
000093,000095,000098,00009F 
000082,000085,000088,00008F 
000093,000095,000098,00009F 
000081,000086,000089,00008E 
000091,000097,00009B,00009C 
000081,000086,00008B,00008C 
000091,000097,00009B,00009C 
000083,000084,00008A,00008D 
000093,000094,00009A,00009D 
000080,000087,000088,00008F 
000092.000095.000099.00009E 


000080,000087 

000090,000097 

000083,000084 

000091,000096 

000082,000085 

000093,000094 

000082,000085 

000092,000095 

000081,000086 

000090,000097 

000081,000086 

000091,000096 

000083,000085 

000093,000094 

000080,000085 

000092,000095 

000081,000087 

000090,000097 

000081,000084 

000091,000096 

000082,000085 

000093,000094 

000082,000085 

000092,000095 

000081,000086 

000090,000097 

000081,000086 

000091,000096 

000083,000085 

000093,000094 

000080,000085 

000092,000095 

000081,000087 

000090,000097 

000081,000084 

000091,000096 


,000089,00008E 
,000099,00009E 
,00008B,00008C 
,00009A,00009D 
,00008A,00008D 
,000098,0Q009F 
,000088,00008F 
,000098,00009F 
,000089,00008E 
,00009B,00009C 
,00008B,00008C 
,00009B,00009C 
,00008A,00008D 
,000098,000090 
,000088,00008F 
,000099,00009F 
,000089,00008E 
,000099,000090 
,00008B,00008C 
,000098,000090 
,00008A,000080 
,000098,000090 
,000088,00008F 
,000099,00009F 
,000089,00008E 
,000099,000090 
,00008B,00008C 
,000098,000090 
,00008A,000080 
,00009A,00009D 
,000088,00008F 
,000099,00009E 
,000089,00008E 
,000099,00009E 
,000088,00008C 
,00009A,000090 



000083,000084,00008A,00008D 
000093,000094,000098,000090 
000080,000087,000088,00008F 
000092,000095,000099,00009F 
000080,000087,000089,00008E 
000090,000097,000099,000090 
000083,000084,OOOC8B,OOO08C 
000091,000096,00009B,00009D 
000083,000085,00008A,00008D 
000093,000094,000098,00009F 
000080,000085,000088,00008F 
000092,000095,000098,00009F 
000081,000087,000089,00008E 
000090,000097,00009B,00009C 
000081,000084,00008B,00008C 
000091,000096,00009B,00009C 


APPENDIX 0. 


V.32 VITERBI DECODER 


;This program is a Viterbi Decoder for V.32. There is a 16 
;time period delay which will approach the maximum possible 
;gain for this type of encoder. If the differential encoder 
;vas not used in the encoder than the call to diff must be 

;commented out. 

• 

f 

{There are 128 memory locations allocated for path memory (8 states x 16 
{time periods sl28 locations). The full scale constellation values are 

{loaded into memory during the intialization routine. 

• 

» 

mfmdecod 

ident 1,1 

page 132,66,3,3,0 

opt nomd,nomex,loc,nocex,mu,cex 

include 'sincos' 
include ’fft’ 


org 

O 

o 

o 

o 



period 

dsm 

128 

{128 locations for path memory 

location 

dsm 

32 

;32 locations for constellation points 

input 

dsm 

16 

;16 locations for input points 

tables 

dsm 

8 

;8 locations for acumulated distance table 

temp 

dsm 

8 

;8 temp locations for distances 

endlong 

equ 

* 


org 

x:endlong 



storr6 

ds 1 



ynow 

ds 4 


;4 locations for input bits 

org 

y;endlong 



ypast 

ds 2 


;2 past bits for differential decoder 


13 boundary tables with 8 points in each of the 4 quadrants 


org 

x:512 


boundryl 

ds 

32 

boundry2 

ds 

32 

boundry3 

ds 

32 

boundry4 

ds 

32 

boundryS 

ds 

32 

boundry6 

ds 

32 

boundry? 

ds 

32 

boundry8 

ds 

32 
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¥ 


boundry9 

ds 

32 


boundrylO 

ds 

32 


boundryll 

ds 

32 


boundryl2 

ds 

32 


boundryl3 

ds 

32 


startfft 

equ 

$100 

;starting location for FFT routine 

points 

equ 

256 

;number of points for FFT 

coef 

equ 

1024 

;location of FFT coefficients 

data 

equ 

1280 

;location of sampled data 

odata 

equ 

1536 

;output data from FFT 

start 

equ 

$250 

;starting location of decoder 

9 

;Load in 

full scale constellation 

values 

9 

four 

equ 

$200000 


three 

equ 

$180000 


two 

equ 

$100000 


one 

equ 

$080000 


zero 

equ 

$000000 


mone 

equ 

$F80000 


mtwo 

equ 

$FOOOOO 


mthree 

equ 

$e80000 


mfour 

equ 

$e00000 


large 

equ 

.9 


small 

equ 

.1 


offset 

equ 

$000000 



sincos points,coef 
opt mex 

org p:startfft 
fft points,data,odata,coef 
opt cex 
org p:start 
jsr initialize 
do #217,_endrun 
jsr readdata 
jsr findmindist 
jsr accumdist 
jsr traceback 
jsr outputdata 
_endrun 

nop 

svi 

t 

; this initialization routine initializes register and 
; modifiers as well as clearing the memeory. The constellation 
; is also loaded into memory here. The accumulated distance 
; array is set so that state zero starts out at a value of 
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; zero and all others start out larger, forcing the paths 

; to merge at the zero states. 

« 
f 

initialize 

move #$ffff,mO 
move #127,ml 
move #$ffff,m2 
move #$ffff,m4 
move #127,m5 
move #15,m6 
move #0,rl 
clr b #$0,r0 
clr a r0,r5 
do #256,clrmem 
move a,x:(r0)+ b,y;(r5)+ 

clriTiim 


move 

#tables+l,r7 

move 
rep #7 

#$400000,al 

move 

al,x:(r7)+ 

move 

#.i npnt ,M 

move 

bl,x:stOir6 

move 

#odata-16,r0 

move 

r0,r2 

move 

#0,x0 

move 

#0,y0 

do 

#16,_clrreg 

move 

x0,x:{r0)+ 

move 

clrreg 

y0,y;(r2)+ 

Nov load 
location. 

full scale values of the constellationin the table 


move 

move 

move 

#location, 

r0,r4 

#mfour,a 

rO 

;Real 

Imag 

move 

move 

#one,b 
a,x: (rO)-f 

b,y:(r4)+ 

; 

1 

move 

move 

move 

#zero,a 
#mthree,b 
a,x:(r0)+ 

b,y:(r4)+ 

: 0 

-3 

move 

move 

#one,b 
a,x: (rO)-f 

b,y:(r4)+ 

; 0 

1 

move 

move 

#four,a 
a,x:(r0)+ 

b,y:{r4)+ 

; 4 

1 

move 

move 

#mone,b 
a,x: (rO)-^ 

b,y:{r4)+ 

; 4 

-1 


;linear addressing 
;modulo 127 addressing 
;linear addressing 
;lii)£i)i‘ addressing 
;modulo 127 addressing 
;modulo 15 addressing 
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f' * 


move 

nzero,a 





move 

//three, b 





move 

a,x:(r0)+ 

b,y:(r4)+ 

> 

0 

3 

move 

//mone, b 





move 

a,x:(r0)+ 

b,y:{r4)+ 

f 

0 

-1 

move 

//mfour,a 





move 

a,x: (r0)-»' 

b,y:(r4)+ 

9 

-4 

-1 

move 

ifmtvo,a 





move 

//three,b 





move 

Ifmona , yl 





move 

a,x: (rO)-f 

b,y:{r4)+ 

9 

-2 

3 

move 

a,x: (roj-f 

yl,y:(r4)+ 

9 

-2 

-1 

move 

//tvo,a 





move 

a,x: (rO)-f 

b,y:(r4)+ 

9 

2 

3 

move 

a,x:{r0)+ 

yl,y:(r4)+ 

9 

2 

-1 

move 

//one, b 





move 

//mthree.yl 





move 

a,x;(r0)+ 

yl,y:(r4)+ 

9 

2 

-3 

move 

a,x:(r0)+ 

b,y:(r4)+ 

9 

2 

1 

move 

//mtvo, a 





move 

a,x:(r0)+ 

yl,y:(r4)+ 

9 

-2 

-3 

move 

a,x: (rO)-f 

b,y:{r4)+ 

> 

-2 

1 

move 

//one, a 





move 

a,xO 





move 

//mthree.a 





move 

//two,b 





move 

b,yO 





move 

//mtvo, b 





move 

a,x:(r0)+ 

b,y:(r4)+ 

9 

-3 

-2 

move 

xO,x:(r0)+ 

b,y:(r4)+ 

9 

1 

-2 

move 

a,x: (rO)'f 

y0,y:{r4) 

9 

-3 

2 

move 

xO,x:(r0)+ 

y0,y;(r4)+ 

9 

1 

2 

move 

//three, a 





move 

a,xO 





move 

//mone, a 





move 

xO,x:(r0)+ 

y0,y:(r4)+ 

9 

3 

2 

move 

a,x:(r0)+ 

y0,y:(r4)+ 

9 

-1 

2 

move 

xO,x;{r0)+ 

b,y:(r4)+ 

9 

3 

-2 

move 

a,x;(r0)+ 

b,y:{r4)+ 

9 

-2 

-1 

move 

//one, a 





move 

//zero,b 





move 

b,yO 





move 

//four,b 





move 

a,x: (rO)-f 

b,y:(r4)+ 

9 

1 

4 

move 

//mthree,xO 





move 

xO,x:{r0)+ 

yO.y:(r4)+ 

9 

-3 

0 

move 

a,x:(r0)+ 

y0,y:(r4)+ 

9 

1 

0 

move 

//mfour,b 





move 

a,x: (rO)-f 

b,y;(r4)+ 

9 

1 

-4 

move 

//mone, a 





move 

a,x:(r0)+ 

b,y:(r4)+ 

9 

-1 

-4 
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mvve 

//three, xO 


move 

x0,x:(r0)+ 

y0,y;{r4)+ 

move 

a,x;(r0)+ 

y0,y:(r4)+ 

move 

//four,b 


move 

a,X t(rO) + 

b,y:(r4)+ 

move 

f!/odata,x0 


move 

//$eff,y0 


move 

x0,y:$424 


move 

rts 

y0,y:$425 



t 

;readdata reads in the data from the output of the FFT. 
i 

readdata 

move y:$/i24,r0 
nop 

move rO,yO 
move //odata-flOO.a 
cmp yO,a 
jseq _delzero 
move x:storr6,r6 
move #>offset,xO 
move x;(rO).a 

add xO,a y:(r0)4,b 

add xO,b a,;<:(r6) 

move b,y(r6) + 

move r6,x:storr6 

move r0,y:$424 

rts 

_delzero 

move #odata+155,r0 
nop 
rts 
i 

;the minimum distance is found to the closest point in every 
; state and stored. The values are stored so that indexing is 
; made easier, state 0,2,3,1,4,7,6,5. This will greatly reduce 
; the number of cycles needed later. A smoothing functi^on is 
; used to accumulate distances in the accumulated table so 

; this minimum distance is multiplied by .1. 

• 

♦ 

;The subroutine findmindist compares the received points to boundaries 
;on the constellation. Once a bounded area is found the closest eight 
;points are read out of the boundary data file and used to update the 
;path distances. 

• 

findmindist 

move x:-(r6),a 
move ^/one,x0 
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cmpm xO,a y:(r6),b 

jgt bigone 

cmpm xO,b //boundryl,r2 

jit continue 

move //tvo,xl 

cmpm xl,b tfboundry4,r2 

jit continue 

move #boundry6,r2 

jmp continue 

bigone move //tvo,xl 
cmpm xl,a 
jgt bigtwo 

cmpm xO,b y/boundi'y2,r2 

jit continue 

cmpm xl,b #boundry5,r2 

jit continue 

bigtwo 

cmpm xO,b #boundry3,r2 

jit continue 

abs a #two,yO 

abs b a,xl 

sub yO,a b,yl 

sub xO,b 

cmpm a,b yl,b 

jgt greateryl 

cmp yO, b boundry?,r2 

jit continue 

move #boundryl2,r2 

jmp continue 
greateryl 

sub yO,b xl,a 

sub xO,a 

cmpm a.b xl,a 

Jgt greatery2 

cmp yO,a l!(boundryl0,r2 

Jit continue 

move yl,b 

cmp yO,b Jl>boundryll,r2 

Jit continue 
move #boundry9,r2 

Jmp centime 
greatery2 

cmp yO,a tfboundryS,r2 

Jit continue 

move #boundryl3,r2 


;x>l 

;x<l,y<l, load r2 with boundry 1 
;and continue 


;x<l,y>landy<2, load r2 with 
;bcundry4, go on 

;x<l,y>2, load r2 with boundry6 
;and continue 


;x>2, Jmp to that case 

;x>l and x<2, y<l load boundryi' 
;and continue 

;x>l,y<2 load boundry 5 and 
;continue 


;x>2 and y<l so load boundry3 
;and continue 
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continue 

clr a x:(r6),xl 

cmp xl,a y:{r6),yl 

jgt negx 

cmp yl,a </24,n2 

jgt posxnegy 
posxposy 

jmp findlst 
posxnegy 

move x:(r2)-fn2,x0 
jmp findist 

negx 

cmp yl,a l!»8,n2 

jgt negxnegy 
negxposy 

move x:(r2)+n2,x0 
jmp findlst 
negxnegy 

move x: (t2)-^n2,x0 
move x: (r2)-fn2,x0 


;update r2 by 24 


; update r2 by 8 


;update r2 by 16 


;The subroutine findist finds the euclidean distance between the received 
;point and the eight points read out of the boundary table. The x and y 
;coordinates are subtracted, squared and added. The square root is not 
;performed. 


iiSt 


move 

x:{r2)+,r0 

move 

//tables, r4 

move 

x:(rO),a 

sub 

xl,a 

sub 

yi.b 

mpy 

xO,xO,a 

mac 

yO,yO,a 

move 

//small, xO 

mpy 

xO,yO,a 

move 

x: (rO),a 

sub 

xl,a 

sub 

yl,b 

mpy 

xO,xO,a 

mac 

yO,yO,a 

move 

//small,xO 

apy 

xO,yO,a 

move 

x:(rO),a 

sub 

xl,a 

sub 

yl,b 

mpy 

xO,xO,a 

mac 

yO,yO,a 

move 

Ismail,xO 

mpy 

xO,yO,a 


y:(rO),b 

a, xO 

b, yO 

x: (r2) + ,r0 
a,yO 

a,y:{r4)+ 

y:(rO),b 

a, xO y:(r4)+,y0 

b, yO 

x: (r2) + ,r0 
a.yO 

y;{r4)+,b 

a,y:(r4)- 

y;(rO),b 

a, xO y:{r4)-,y0 

b. yO 

x:(r2>+,r0 
a,yO 
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move x:(rO) ,a 
sub xl,a 
sub yl,b 
mpy xO,xO,a 
mac yO,yO,a 
move //small,xO 
mpy xO,yO,a 
move x:(rO),a 
sub xl,a 
sub yl,b 
mpy xO,xO,a 
mac yO,yO,a 
move jir'small,xO 
mpy xO,yO,a 
move x:(rO),a 
sub xl,a 
sub yl,b 
mpy xO,xO,a 
mac yO,yO,a 
move //small,xO 
mpy xO,yO,a 
move x:(rO),a 
sub xl,a 
sub yl,b 
mpy xO,xO,a 
mac yO,yO,a 
move //small,xO 
mpy xO,yO,a 
move x:(rO),a 
sub xl,a 
sub yl,b 
mpy xO,xO,a 
mac yO,yO,a 
move //small,xO 
mpy xO,yO,a 
move a,y:(r4) 
rts 


a,y:(r4)+ 
y:(rO),b 

a, xO 

b, yO 

x:(r2)+,r0 
a,yO 

a,y;(r4)+ 

y:(rO),b 

a, xO y:{r4)4,y0 

b, yO 

x: (r2) + ,r0 
a,yO 

a,y:(r4)+ 
y;(rO),b 

a, xO y:(r4)+,y0 

b, yO 

x;(r2)+,r0 
a,yO 

y: (r4) + ,b 
a,y:(r4)- 
y;(rO),b 

a, xO 

b, yO 

x:(r2)+,r0 
a,yO 

a,y:(r4)- 

y:(rO),b 

a, xO 

b, yO 

x:(r2)+,r0 
a,yO 


;the accumulted distance routine adds the smallest distance 
;from the previously computed table for all pathes going 
;a state and does this for all eight states. 

accuffldist 


clr a 

//tables,rO 

move 

//$7fffff,al 

move 

rO, r4 

move 

Hemp , r2 

move 

//3,m0 

move 

m0,m4 

move 

(!/2,nl 

move 

nl,n5 
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move 


rl,r5 


;Distances in the accumulated distance table are added to distances in the 
;path table and compared for the four paths. This is done by incrementing 
;through a specially ordered path table. 

;find minimum distance to state zero 
do )/4,statezero 

move x:(rO),xO y:(r4),b 

add xO,b 

cmp b,a 

tge b,a r0,r3 

tge b,a r4,r7 

move x:(rO)+,xO y;(r4)+,b 

statezero 

move r3,x:(rl)+nl 
move a,x:{r2)+ y:(r4)+,b 

clr a r7,y:{r5)+n5 

move il>$7fffff,al 


;find minimum distance to state two 
do #4,statetwo 
move x:(rO),xO y:(r4),b 

add xO,b 
cmp b,a 

tge b,a r0,r3 

tge b,a r4,r7 

move x;{rO)+,xO y:(r4)-,b 
statetvo 

move r3,x:(rl)+ni 
move a,x:(r2)+ y:(r4)+,b 

clr a r7,y:(r5)+n5 

move i'/$7fffff,al 

;find minimum distance to state four 
do {!^4,statefour 
move x:(rO),x0 y:(r4),b 

add xO,b 
cmp b, a 

tge b,a r0,r3 

tge b,a r4,r7 

move x:(rO)+,xO y:(r4)4,b 

statefour 

move r3,x;(rl)4nl 
move a,x:(r2)4 y:(r4)4,b 

clr a r7,y:(r5)4n5 

move )!>$7fffff,al 

;find minimum distance to state six 
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do /^4,statezsix 
move x;{rO),xO y:(r4),b 

add xO,b 
cmp b,a 

tge b,a r0,r3 

tge b«a r4,r7 

move x;(rO)+,xO y:{r4)-,b 

statezsix 

move r3,x:(rl)-nl 
move atx:(r2)-f 
move r7,y:(r5) 
move i^tables-»'4,r4 
move r4,r0 
move x:(rl)-nl,a 
clr a x:(rl)-,b 
move //$7fffff,al 
move rl,r5 

;find minimum distance to state one 
do //4,stateone 
move X!(rO),x0 y:(r4),b 

add xO,b 
cmp b,a 

tge b,a r0,r3 
tge b,a r4,r7 
move x;(rO)+,xO y:(r4)+,b 

stateone 

move r3,x:(rl)+nl 
move a,xj(r2)+ y:(r4)+,b 

clr a r7,y:(r5)+n5 
move #$7fffff,al 

;find minimum distance to state three 
do f/4,statethree 
move x:(rO),x0 y:(r4),b 

add xO,b 
cmp b,a 

tge b,a r0,r3 
tge b,a r4,r7 
move X:(rO)+,xO y:(r4)-,b 

statethree 

move r3,x; (rl)4^nl 
move a,x:(r2)-f y:(r4) + ,b 

clr a r7,y;(r5)+n5 
move #$7fffff,al 
move (r4)+ 

;find minimum distance to state five 
do ll^4,$tatefive 
move x:(rO),xO y:(r4),b 

add xO,b 
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cmp b,a 

tge b,a rO,r3 

tge b,a r4,r7 

move x:(rO)+,xO y:{r4)-,b 
statefive 

move r3,x:(rl)+nl 
move a,x:(r2)-t- y:(r4)-,b 

clr a r7,y:{r5)+n5 
move tf$7fffff,al 

;find minimum distance to state seven 
do #4,stateseven 
move x:(rO),x0 y:(r4),b 

add xO,b 
cmp b,a 

tge b,a r0,r3 
tge b,a r4,r7 
move x:(rO)-<-,xO y:(r4) + ,b 

stateseven 

move r3,x:(rl)+ 
move a,x:(r2)+ y:(r4)+,b 

clr b r7»y:{r5)+ 
move #$7fffff,bl 

;nov move nev accumulated distances into the accumulated 
;distance table from the temporary table also find the min 
;distance state and store in r4 which is no longer used 

move i!»$ffff,mO 
move //$ffff,m4 
move <^temp,r3 
move )i>tables,rO 
move il»large,xl 
move #2,n0 
do i!>4 ,endtable 
move x:(r3)+,x0 
mpy xl,xO,a 

cmp a,b a,x:(rO)-«-nO 

tge a,b r0,r4 

endtable 

move #tables-^l,rO 
do #4,endtablex 
move x:(r3)+,x0 
mpy xl,xO,a 

cmp a,b a,x:(rOj-fnO 

tge a,b r0,r4 

endtablex 
> 

:store in rO instead of r4 
i 

move r4, rO 


80 





move 

move 

rts 


//8,nl 

(rO)-nO 


;the traceback routine now goes back through every time period 
;8tarting with the current time period and finds the state 
;from which the path came from one time period previous. At 
;the end of this search, the last state found will also point 
;to the path at that state, which is the output of the 
;trellis. 

traceback 

;find the displacement from the pointer to table and store 
;value in n4 

move f/tables,nO 
move {rl)-nl 
lua (rO)-nO,n5 
move rl,r5 
do //15,endtrace 
move {rl)-nl 
move x:(r5+n5),r0 
move rl,r5 
lua (r0)-n0,n5 
endtrace 

move //location,rO 
move y:(r5+n5),a 
rts 

;the output data routine unscrambles the path order and finds 
;one of the four points on the constellation coresponding to 
;the output state which is closest to the original input at 
;that time period. 

outputdata 


move 

a,b 



move 


//>$bl 

,x0 

cmp 

xO,a 

#>$b2. 

yo 

teq 

yo,b 



cmp 

yO,a 

</>$b3. 

xO 

teq 

xO,b 



cmp 

xO,a 

#>$bl. 

yo 

teq 

yO,b 



move 


#>$b5, 

,x0 

cmp 

xO,a 

//>$b7, 

yo 

teq 

yo,b 



cmp 

yO,a 



teq 

xO,b 



move 

!ii),r2 
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¥ 


move 

^/tables,n2 

move 

x:storr6,r6 

lua 

(r2)-n2,n3 

move 

n3,a 

asl 

a 

asl 

a 

move 

a,nO 

move 

r6,r3 

lua 

(r0)+n0,r4 

move 

#>$7fffff,xl 

move 

r4,r0 

do 

f!f 4, endout 

move 

x: (r3),a 

move 

x: (rO) + ,xO 

sub 

xO,a 

sub 

yO,b a,xO 

mpy 

xO,xO,a b,yO 

mac 

yO.yO.a 

tfr 

a,b xl,a 

cmp 

xl ,b 

tit 

b,a r0,r7 

move 

a,xl 

endout 


clr 

a (r7)- 

move 

f/location.nO 

move 

r7,r0 

move 

//$f,al 

lua 

{r0)-n0,r7 

move 

r7,x0 

and 

xO,a 

jsr 

diff 

move 

y:$425,r0 

nop 


move 

aO,y:(r0)+ 


move r0,y:$425 
rts 


y:(r6),b 

y:(r4)+,y0 


;The subroutine diff differentially decodes the two most significant bits. 
;Each bit is stored in its own memory and the bits are decoded using 
;Qln = Yin EOR Yln-1, Q2n = (Qln AND Yln-1) EOR Y2n-1 EOR Y2n. The four 
:bit symbol is formed and output. 

t 

kff 

move //ynow+3,r0 
move 

move #ypast,r7 
do Ilf 4, diff loop 1 
and xO,a a,xl 
move al»x:(rO)- 

move xl,a 
asr a 
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diffloopl 


move 

x:(rO)+,a y:(r7)+,y0 

move 

x:{rO)+,a y:(r7)-,yl 

move 

x:(rO)-,b a,y:(r7)+ 

move 

b,y:(r7)- 

eor 

yO,a a,xO 

eor 

yl,b a,xl 

and 

yO.a b,yl 

eor 

yl,a xl,b 

move 

b,x:(r0)+ 

move 

a,x:(rO)- 

clr 

a 

clr 

b 

do 

#4,diff2 

move 

x:(r0)+,b0 

addl 

b,a 


diff2 

rts 
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